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ABSTRACT 

Context. Based on color-selected IRAS point sources, we have started to conduct a survey of 47 high-mass star-forming regions in 
the southern hemisphere in 870 fim dust continuum and molecular line emission in several frequency ranges between 290 GHz and 
806 GHz. This paper describes the pilot study of the three sources IRAS 12326-6245, IRAS 16060-5 146, and IRAS 16065-5 158. 
Aims. To characterize the physical and chemical properties of southern massive star-forming regions, three high-luminosity southern 
hemisphere hot cores are observed in continuum and molecular line emission. Based on the results obtained in the three sources, 
which served as templates for the survey, the most promising (and feasible) frequency setups for the remaining 44 sources were 
decided upon. 

Methods. The sources were observed with the Atacama Pathfinder Experiment (APEX) in five frequency setups aimed at groups 
of lines from the following molecules: CH 3 OH, H 2 CO, and CH 3 CN. Using the LTE approximation, temperatures, source sizes, and 
column densities were determined through modeling of synthetic spectra with the XCLASS program. Dust continuum observations 
were done with the Large APEX BOlometer CAmera (LABOCA) at 870 //m and the 3 mm continuum was imaged with the Australia 
Telescope Compact Array (ATCA). 

Results. Based on the detection of high-excitation CH 3 CN lines and lines from complex organic species, the three sources are 
classified as line rich, hot core type sources. For all three, the modeling indicates that the line emission emerges from a combination 
of an extended, cooler envelope, and a hot compact component. All three sources show an overabundance of oxygen-bearing species 
compared to nitrogen-bearing species. While the chemistry in the three sources indicates that they are already in an evolved stage, 
the non-detection of infrared heating sources at the dust continuum peak of IRAS 16065-5 158 points to this source still being deeply 
embedded. Because this work served as a pilot study, the approach to observe the remaining 44 massive star-forming regions was 
chosen based on its results. 

Conclusions. The three sources are massive, luminous hot cores. While IRAS 16065-5 158 seems to be a very young deeply embedded 
object, IRAS 12326-6245 and IRAS 16060-5 146 seem more evolved and have already developed UCH11 regions. 
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1. Introduction 

Massive stars (M > 8 M Q ) are important building blocks of 
galaxies. Although they are rare, they provide most of the stellar 
luminosity and heavily influence their surroundings through en- 
ergetic interactions throughout their short lifetime. A further un- 
derstanding of the role of massive stars involves understanding 
their formation processes, which has always posed a challenge, 
because the earliest stages of massive star formation are deeply 
embedded in these objects' natal dust clouds. 

One approach to study high-mass star-forming regions in 
depth would be to observe a large number of sources in a wide 
frequency range, which is not yet feasible. Observing the sources 
in a number of frequency setups to cover many molecular lines 
is needed to study the prevailing excitation conditions of the 
molecular gas. In the past, individual sou rces have either been 
studied in several frequency setups (see ICesaroni et al 



1994 



2006, 



iMacDonald eTaT] 119961: lOlmi et al.l |2003; Beuther etal 
to name a few), or surveys of several sources have been per- 
formed, mostly on the no rthern hemisphere, with only very lim- 
ited f r equency coverage (|Churchwell et al.l l 1992c iMolinari et al 



200 1 iMuelleretaD l2002t ISridharan et alJ l2002t IShirlev et al 
2003l) . On the southern hemisphere, iPurcell et al. (120061) and 



lUrquhart et al.l (l2007bl) have performed line surveys in the 
mm regime. While many samples of high-mass star-forming 
regions have been observed in the dust continuum on the 
northern and southe r n hemispheres , such a s the sam p les of 
iBeuther et all d2002l): iFaundez et all d2004l) : iHill et all d2005l) 
and lBeltran et al.l (120061) . no project has so far observed the same 
sample both in a high-frequency line survey and in dust contin- 



To address the lack of statistically significant sub-mm sur- 
veys of southern hemisphere sources in both lines and contin- 
uum, we have conducted a survey of 47 high-mass star-forming 
regions in the southern hemisphere in 870 fim dust continuum 
and molecular line emission in several frequency ranges between 
290 GHz and 806 GHz, based on color-selected IRAS point 
source criteria. Our se lection criteria follow those described in 
ISridharan et alJ d2002l) . with the exception that our sample in- 
cludes sources with and without detected radio continuum emis- 
sion in medium resolution single dish 5 GHz surveys. 
This large sample contains sources at very different stages of 
evolution and will allow a consistent and comparative analysis 
of these regions. 
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Out of this survey, we will present in this paper observa- 
tions of three high-luminosity southern hemisphere hot cores 
with the Atacama Pathfinder Experiment telescope, which 
served as targets for a template study for our large survey. 
IRAS 12326-6245, IRAS 16060-5146, and IRAS 16065-5158 
(TableQ] from now on we omit the IRAS in the name for brevity) 
were chosen to be analyzed in this pilot study, because they 
are the most luminous and line-rich sources in the sample of 
our study based on their HCO + and SO2 spectra, and therefore 
the most suitable targets to test the feasibility of the project. 
A fourth very strong hot molecular core of the sample, IRAS 
17233-3606, has been analyzed and the results have been pub- 
lished separately dLeurini et al.ll2008T) . Based on the results ob- 
tained in these three sources, the most promising (and feasible) 
frequency setups for the remaining 44 sources were selected. 



1.1. Source details 

The three sources had alr eady been observed in the course 
of several maser st udies JWalsh et all JlWfr ICaswelll Il998t 
iMacLeod et aIlll998HWalsh et al.ll 1 998h ICaswellll20q lh. as well 
as in the m m continuum survey of IFaundez et al.l (12004b . the 
CS survey of iBronfman et al.l (1 19961) and the HNCO survey of 
IZinchenko et al.l d2000h ~ 

We derived the kinematic distances to the sources from the 
center velocities of the C n O(3-2) line by usin g the Galactic 
velocity field description of lKalberla et all (120071) . In the follow- 
ing paragraph, the characteristics of the sources are summarized. 



IRAS 12326-6245 

The kinematic distance estimated for 12326-6245 from its 
systemic velocity of -39.3 kms -1 is 4.4 kpc. 12326-6245 
features cl ass II CH3OH masers at 6.7 GHz, OH masers at 
1.6 GHz (ICaswelll 1 19981) and H 2 masers (IMacLeod et alJ 
[19981 Observed at 1 .2 mm, 12326-6245 shows compact dust 
continuum emission dHenning et al.| | 2000t IFaundez et al.l 120041 
iHill et al.ll2005b IMiettinen et al.ll2006l) . which coinci des with the 
location of two peaks of cm continuum emission ((Walsh et al. 
119981: lUrquhart et al.|[2007aT) observed at high resolution. These 
peaks correspond to the location of two strong MIR peaks a t 
the Q (20 fim) and N (10 fim) bands dHenninget al.l 120001) . 
Molecular line observations in the CH 3 CN(5^fl CH 3 CN(6- 
5), CH 3 CN(8-7) and CH 3 CN(12-11) series reveal a two- 
component temperature st ructure hinting at a hot core embed- 
ded in a cooler envelope (lArava et al.l 12005b . Maps of molecu- 
lar line data show that 12326-6245 seems to have one of the 
largest bipolar outflows detected in th e southern hem isphere 
dHenning et al.l 120001: IWu et al.l 120041) . iMiettinen et ail d2006l) 
derive a rotation temperature of 34 K from observations of 
CH 3 CCH bands. 



1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory and the Onsala Space Observatory. 

2 Here and below we designate series of lines close together in fre- 
quency. The first and second numbers are upper and lower values of 
the total angular momentum quantum number J of the individual transi- 
tions of the series that have different projected angular momentum (A*) 
values. 



IRAS 16060-5146 

The source has class II CH 3 OH mase rs at 6.7 GHz and 
OH masers at 1665 MH z, 1667 MHz (ICaswelll [19981) and 
6035 MHz (ICaswellll200l[). It shows strong, spherical 1.2 mm 
dust continuum emission (IFaundez et alJl2004t IHill etal]l2005t 
IMiettinen et al.l 120061) that peaks at the location of the maser 
sources. H igh-resolution radio c ontinuum observations at 3.6 
and 6 cm (lUrquhart et al.l l2007al) show two UCHn regions at 
the location of the mm peak. We derived 5.5 kpc and 9.4 kpc 
for the near and far distances, respectively, from the systemic 
velocity of 16060-5146 (-90.95 km s -1 ) found in our data. The 
kinematic distance estimates for this source found in the litera- 
ture differ somewhat between 5.3-6.1 kpc for the near distance 
and 9.6 kpc for the far distance, which can be explained with the 
spread found in the systemic velocities assumed for this source, 
which range from -88 km s to -91 kms -1 . 

IRAS 16065-5158 

In 16065-5158, OH maser emission at 1665 MHz, 1667 MHz, 
and 16 12 MHz and a CH 3 OH 6.7 GHz maser has been ob- 
served dCaswelll [l998h. 16065-5158 sho ws extended 1.2 mm 
continuum emission (IFaundez et al.ll2004l) and a UCHn region 
dUrquhart et al.ll2007al) . We derive a kinematic distance of 4 kpc 
and 10.8 kpc for the near and far values of 16065-5158, using a 
systemic velocity of -62 . 1 7 km s~ 1 . 



2. Observations and data reduction 

2.1. Atacama Pathfinder Experiment 
2.1.1. Line data 

The sub-mm data were taken between November 2005 and 
November 2 006 wi th th e Atacama Pathfinder Experiment 
iGiisten et al.l (APEX. |2006|) telescope on Chajna ntor in Chile. 
We u sed both the APEX 2a heterodyne receiver dRisacher et al.l 
I2006I) and the dual frequency 460/810 GHz Fi rst Light APEX 
Subm illimeter Heterodyne (FLASH) receiver dHevminck et al.l 
120061) for our observations (see Table|2]for details on the receiver 
setups). As backend, the MPIfR-built F ast Fourier Transform 
Spectrometer (FFTS) dKlein et al. 1 120061) was used, which con- 
sists of two units with 1 GHz band width each. The pointin g ac- 
curacy was checked on G327. 3-0.6 dWvrowski et al.ll2006l) . and 
was found to be accurate within 2'.'. Most of the data were ob- 
served as single pointings. When mapping was done, raster maps 
with a map size of 407 x 407 and beam spacing were observed. 
The observations were done in the position-switch mode. The 
off-positions were chosen with an offset of (600'.', 07) in azimuth 
and elevation from the pointing position and were checked to be 
free of CO emission. All observations were obtained in double 
sideband mode. Calibration errors were estimated to be on the 
order of 30%. These errors are influenced by pointing and atmo- 
spheric fluctuations, but the largest part of the uncertainty stems 
from uncertain knowledge of the sideband gain ratios. The data 
were reduced with the CLASS software package- 
Table |2] lists the frequencies of the bands of the molecular 
lines observed with APEX that are subject of the pilot study. 
These lines were the targets of the frequency setups, but because 
up to 1.8 GHz bandwidth could be used, lines from many other 
molecules could be observed simultaneously. 



http://www.iram.fr/IRAMFR/GILDAS 
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Table 1. Peak coordinates of the 1.2 mm dust continuum data, distances and IRAS luminosities of the three hot core sources (see 
text). 



Source R.A. Dec. l,b d (near/far) L (near/far) 

(J200Q) (J2000) (kpc) (1O 5 L ) 

IRAS 12326-6245 12:35:35.90 -63:02:29.00 301.13783,-0.22483 4.4 2.7 

IRAS 16065-5158 16:10:20.01 -52:06:13.26 330.87748,-0.36821 4.0/ 10.8 3.0/21.2 

IRAS 16060-5146 16:09:51.40 -51:55:06.98 330.94914,-0.18246 5.5/9.4 7.8/23.8 



Table 2. Frequency setup for the molecular line observations. 



Transition 


Center Frequency 


Receiver 


©HPBW 




B e ff 


Tracer 


T 

1 sys 




(GHz) 




('■') 








(K) 


CH 3 OH(6-5) 


290.8 


APEX2a 


20 


0.95 


0.73 


T 


170- 190 


CH 3 OH(7-6) 


337.8 


APEX2a 


17.5 


0.95 


0.73 


T 


180-210 


CH 3 CN(16-15) 


294.5 


APEX2a 


20 


0.95 


0.73 


T 


140 - 200 


H 2 CO(4-3) 


290.0 


APEX2a 


20 


0.95 


0.73 


T 


214-291 


H 2 CO(6-5) 


436.9 


FLASH460 


14 


0.95 


0.6 


T 


730- 839 


CO (3-2) 


345.8 


APEX2a 


17 


0.95 


0.73 





550 


CO (4-3) 


461.0 


FLASH460 


13 


0.95 


0.6 





1000-1100 


CO (7-6) 


806.7 


FLASH810 


7 


0.95 


0.43 





5100-5700 


C n O(3-2) 


337.8 


APEX2a 


17.5 


0.95 


0.73 


cd 


180-210 


HCO + (4-3) 


357.0 


APEX2a 


17 


0.95 


0.73 





240 - 324 



Notes. In the column marked 'Tracer', T stands for temperature, o for outflow and cd for molecular column density. 



To analyze the chemical composition and physical properties 
of the three sources, five major setups were used, CH 3 OH(6- 
5) at 291 GHz (together with H 2 CO(4-3) for 16060-5146 and 
16065-5158), CH 3 OH(7-6) at 338 GHz, CH 3 CN(16-15) at 
294 GHz, H 2 CO(4-3) at 291 GHz and H 2 CO(6-5) at 436 GHz. 
The positions of the hot cores were first defined as the peak posi- 
tions of HCO + (4-3) maps at 357 GHz and later refined by APEX 
continuum scans. 

On-source integration times for APEX2a were around 5-8 min 
per single pointing setup, and 12-20 min per pointing for 
FLASH observations, depending on the weather conditions. 

2.1.2. Continuum data 

In summer 2007, we obtained continuum data for our sources 
from the Large Apex BOlometer Camera (LABOCA). 
LABOCA is a 295 channel bolometer at 87 yum with a field o f 
view of 1 1 .4 ' and a beam FWHM of 18.6 '/ (ISiringo et alj2009h . 
The pixel separation is 36 To produce a fully sampled map, the 
data were observed in the raster-spiral mode, with 35 s integra- 
tion time on source, a spacing of 27 '.' and 4 subscans. This sets 
the scanning velocity to about 4'/s and the total integration time 
on source was on average 2.3 min. 12326-6245 was observed in 
two scans. 

The atmospheric opacity was determined every hour through 

skydips. The zenith opacity was 0.3 for 12326-6245 and 

0.09 for 16060-5146 and 16065-5158. The rms noise for 

12326-6245 is 100 mJy/beam, while the rms noise for 

16060-5146 and 16065-5158 is 50 mJy/beam. 

IRAS 1 6293-2422, IRAS 1 3 1 34-6242 and i] Carinae were used 

as flux calibrators and the pointing was checked on 77 Carinae 

for 12326-6245 and on IRAS 16293-2422 for 16060-5 146 and 

16065-5158. The pointing was good within 3'.'. 

The data were reduced with the BoA software package (Schuller 

et al., in prep.). 



2.2. ATCA 

High spatial resolution data at 3 mm were taken in October 2006 
with the Australia Telescope Compact Array (ATCA) interfer- 
ometer at Narrabri, Australia. The array was in the H75 con- 
figuration with a synthesized beam of &.' x 4'.' for 16060-5146, 
5'.' x 5'.' for 16065-5158 and &.' x 6'.' for 12326-6245 and in the 
FULL_16_256-128 correlator configuration. 
The continuum was observed at 87.89 GHz with a bandwidth of 
128 MHz and a resolution of 8.8 MHz. 

The sources were observed during two nights, PKS 1253-055 
was used as bandpass calibrator, Mars was observed on 
the second night as flux calibrator, the gain calibration for 
12326-6245 was done on PKS 1147-6753, and on PKS 
1613-586 for 16060-5146 and 16065-5158. The theoretical 
continuum rms noise level was 3 mJy/beam for the continuum 
and 130 mJy/beam for the line observations, respectively. The 
data were reduced with MIRIAE0 and t he dirty images were de- 
convolved with the CLEAN algorithm dHogbomll974l) . 

3. Observational results 

3.1. Line data 

Because in the beginning LABOCA data were not available, the 
three sources were mapped in HCO + (4-3) at 357 GHz to ver- 
ify A6_P^^rJO^tion_observed in the 1.2 mm continuum maps 
of iFaundez et al.l d2004h . First, the regions were mapped with 
beam spacing to locate the peaks, then the coverage was refined 
with fully sampled maps around the peak positions (see Fig. [T}. 
The peak positions in 12326-6245 and 16060-5146 are signif- 
icantly offse t from the 1.2 mm co ntinuum peak positions as de- 
termined bv lFaundez et al. (2004). This is mostly likely owing 
to a pointing problem in the 1.2 mm continuum data, because in 
12326-6245 and 16060-5146, respectively, the peak positions 

4 http://www.atnf.csiro.au/computing/software/miriad 
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derived by iFaundez et al.l (12004]) are offset by 6 and 15 '.' from 
the LABOCA and ATCA peak positions, which agree well for 
all three sources. The final pointing positions of the line obser- 
vations can be seen in Table [3] The offset relative to the ATCA 
positions is less than a FWHM in all cases. When comparing re- 
sults of different frequency setups later on, note that the offsets 
between the different positions are less than 1/3 of the FWHM 
in all but one cases. The extent of the sources in HCO + (4-3) is 
larger than the beam, which justifies the determination of beam- 
averaged column density later on. 

In order to determine the total column density and the prop- 
erties of the envelopes, the sources were also mapped in CO(3- 
2) (Fig. ID and observed in CO(4-3), CO(7-6) and C n O(3- 
2) on the peak positions (see Fig. |3J. Both the CO lines and 
the HCO + (4-3) line show clear signs of outflow activity in 
the wings of all sources and self-absorption at the velocity 
of C I7 0. All three sources have very high-velocity outflows 
with a width of the blue wing, Avbi ue , of 30km s and a 
red wing width, Av le d, of 20km s in both 16060-5146 and 
16065-5158. In 12326-6245, which was alread y described as 
one of the most massive and energetic outflows bv lHenning et alj 
d2000l) . we found Ai'biue = 50km s _1 and Av re( j = 60km s -1 . 
16060-5146 and 16065-5158 show an asymmetric line pro- 
file with a stronger blue shoulder at the peak position of the 
HCO + (4-3) emission, which is a sign of motion in the core. At 
the resolution of the current study, it is not possible to distin- 
guish between infall and rotation as possible origin of the large 
scale motions. Note that in both 12326-6245 and 16065-5158, 
the high-/ CO the intensity in the outflow wings is stronger than 
in the lower-/ lines and one sees stronger self-absorption. In 
16060-5146, however, the outflow wings are less strong in the 
CO(7-6) line, indication of a lower CO abundance in the high- 
velocity gas. 

The three sources were observed in CH30H(6-5), 
CH 3 OH(7-6), CH 3 CN(16-15), H 2 CO(4-3), and H 2 CO(6- 
5). As will be further discussed in the next section, all 
three sources exhibit a rich spectrum of molecular lines. In 
12326-6245 and 16065-5158, we found 18 different species, 
while in 16060-5146, we found 16 species. As is typical for hot 
core sources, high-excitation lines up to energies of 346 K above 
the ground state and lines from vibrationally excited states were 
observed. Especially the setup around 338 GHz is very rich 
in lines, with detections of up to 40 lines per G Hz. The lines 
were identified with the programme XCLASS (ISchilke et alJ 
I1999L IComito et al.ll2005h and the rest frequencies provided by 
the Cologne database for Molecular Spectroscopy (CDMSjj 
dMiilleretalJ l200lL 120051) a nd JPL molecular spectroscopy 
databas d3 ( Pickett etalJll998l) . Because the line identification 
went hand in hand with the LTE modeling of the lines, it is 
described in more detail in Sec. 14.11 To give an example of 
the chemical complexity of the sources, the identified lines 
that were included in the XCLASS modeling are listed for the 
source 16065-5158 in the online material. Note that it was not 
possible to perform individual Gaussian fits to all the CH3OH 
and CH3CN lines included in this list, because some of them 
were heavily blended with lines from other species in all three 
sources. 

For CH3OH, CH3CN, and H 2 CO, the parameters of the lines 
(determined from Gaussian fits) can be found in Tablesl4ltoll2lin 
the online material. The average line widths of the CH3OH lines 
in the three sources are 6 km s~ , 10 km s~ and 7 km s~ and of 
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Fig.l. Maps of the HCO + (4-3) line for 12326-6245 (top), 
16060-5146 (center) and 16065-5158 (bottom). The (0,0) posi- 
tions correspond to the 1 .2 mm peak positions of Faundez (2004) 
listed in Tab.Q] At each offset position, the velocity and T m b tem- 
perature scales go from -100 - 20 kms~' and -5 - 25 K, -150 - 
km s" 1 and -5 - 20 K and -120 - km s" 1 and -5 - 20 K re- 
spectively for 12326-6245, 16060-5146 and 16065-5158. 



the CH 3 CN(16-15) lines 6 km s _1 , 11 km s" 1 and 8 km s" 1 for 
12326-6245, 16060-5146 and 16065-5158 respectively. 

16060-5146 has a slightly broader profile than the other 
two sources. In both CH?CN(16-15) and the high-excitation 
CH30H(v,=l) lines in 16060-5146, one can see a double- 
peaked profile (Fig. [4] and Fig. [T3]), which might suggest that 
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Table 3. Coordinates of the line observations of the three sources. 



Source 




R.A. 


Dec. 


Frequency Setup 






\JZ\J\J\J) 


f'TinniYi 

{JZ.\J\J\J) 




IRAS 12326- 


-6245 


12:35:33.7 


-63:02:19.0 


H 2 CO(4-3) 


IRAS 12326- 


-6245 


12:35:34.4 


-63:02:29.0 


CH 3 CN(16-15),CH 3 OH(6-5) 


IRAS 12326- 


-6245 


12:35:35.9 


-63:02:29.0 


CH 3 OH(7-6) 


IRAS 16065- 


-5158 


16:10:20.0 


-52:06:03.3 


CH 3 CN(16-15),CH 3 OH(6-5) 


IRAS 16065- 


-5158 


16:10:20.6 


-52:06:08.3 


H 2 CO(4-3), H 2 CO(6-5), CH 3 OH(7-6), high-/ CO 


IRAS 16060- 


-5146 


16:09:53.6 


-51:55:07.0 


CH 3 CN(16-15),CH 3 OH(6-5) 


IRAS 16060- 


-5146 


16:09:53.9 


-51:54:59.5 


H 2 CO(4-3), H 2 CO(6-5), CH 3 OH(7-6), high-/ CO 



the larger line widths seen in this source in the other species are 
caused by the superposition of two velocity components. 

In 16065-5158, a similar double-peaked profile can be 
found for the highly (torsionally) excited CH 3 OH(v,= l) lines. 
16065-5158 also displays strong CH 3 CN lines from the vs = 1 
bending mode, which is not seen toward the other two sources. 
Looking at the distribution of peak velocities versus energy of 
the lines, one can see no trend in all three sources and molecules. 
The peak velocities of 16060-5146 are scattered over a broader 
field than those of the other two sources. 

3.2. Continuum data 
3.2.1. LABOCA 

While 12326-6245 (see Fig. [5] for the three sources) shows a 
spherical, centrally peaked morphology at 870 /im, 16060-5146 
and 16065-5158 show an extended morphology with one and 
three peaks, respectively. The latter two sources are also sur- 
rounded by secondary peaks, unlike 12326-6245, which is 
isolated within the LABOCA field of view. Here, the bright 
MIR cluster seen in the Spitzer Space Telescope Galactic 
Legacy Infrared Midp lane Survey Extraordinaire (GLIMPSE) 
(Benjami n et alj|2003l) emission, on which 12326-6245 is cen- 
tered is the only sign of activity in the region, while the re- 
gion immediately around the dust emission traced by LABOCA 
is dominated by enhanced 8 fim emission, a tracer of photon- 
dominated regions (PDRs). 

In 16060-5146, the dominant feature in the LABOCA field 
of view is a large bubble that is visible in the MIR emission. 
The 870 fim emission of 16060-5146 corresponds to an active 
star-forming site located, at least in projection, at the edge of the 
bubble. The other cores traced by 870 fim emission can be found 
tracing the rim of the bubble and in another region to the south 
of 16060-5146. 

16065-5158, which is only 12' away from 16060-5146, is 
located in an environment with very active star formation. One 
finds a very extended region of 8 fim emission, interspersed with 
infrared dark clouds and an association of bright young stars at 
the center. 16065-5158 is located right at the center of the asso- 
ciation, in a deeply embedded region free of MIR emission. An 
elongated structure at 4.5 fim (color-coded green) seems to be 
coming from this deeply embedd ed region. Because 4.5 urn is 
associated with hot shocked gas dCvganowski et al.1 l2008). this 
situation is very suggestive of a massive outflow stemming from 
16065-5158. The very extended 870 fim emission found in the 
LABOCA field of view for 16065-5158 traces nearly the whole 
extent of the active star-forming region visible in the GLIMPSE 
data. 



Table 13. Properties of the 870 fim continuum as derived from 
the LABOCA observations. 



Source 




S peak 


S int 


AT(H 2 ) 






(Jy/beam) 


(Jy) 


(xl0 23 cirT 2 ) 


12326- 


-6245 


19.9 


42.6 


4.2 


16065- 


■5158 


17.9 


41.8 


3.8 


16060- 


-5146 


46.2 


81.0 


8.2 



To obtain source positions, sizes, and peak fluxes, the 
LABOCA maps (see Fig.[5]l were analyzed with the sfind routine 
in the MIRIAD software package. The integrated fluxes were ob- 
tained over the area inside the contour line representing 10% of 
the peak flux, which corresponds to an average signal-to-noise 
ratio of 120. The size of the sources was determined fitting a 
circularly symmetric two dimensional Gaussian source model. 
The source positions, fluxes, sizes, and the calculated H 2 column 
density can be found in Tables [13] and [14] The positions agree 
within the pointing uncertainties with the 3 mm peak positions 
derived with ATCA. The column density was ob t ained from Eq. 
[TJ following lMotte et alJ (ll998h : lMiettinen et al.1 d2006l) 



rdust 



JV(H 2 ) = 1.67 x 10 



-22 



ByiT^fmiHKdRa 



(1) 



where I v is the peak flux in Jy/beam, fi v (Td) the Planck func- 
tion of a blackbody at dust temperature r<j , fi is the mean molec- 
ular mass assuming 10% contribution of helium, kj the dust 
absorption coefficient of 0.176 m 2 /kg linearly approximated at 
87 iim from model V (th i n ice m antles, n = 10 6 crrr 3 ,/3 = 1.8) 
of lOssenkopf & Hennin g| ill 994 and a dust-to-gas mass ratio 
R& = :jir. Model V is suited for sources in which consider- 
able depletion on ice occurs. Unlike the models with thick ice 
mantles, which apply for dark clouds, the conditions in model V 
include the in fluence of a heating sou r ce tha t is starting to evap- 
orate the ices. lOssenkopf & Hennind (fT994) give the uncertain- 
ties for the dust absorption coefficient to be a factor of 2 for ice 
covered dust and state that it can be up to a factor of 5 higher in 
disk regions, where the ice mantles are already evaporated off the 
grains. As provided by model V, we use ft — 1 .8 (k ~ (v/vof) in 
this work, which is consistent w ith values of (3 = 1 .5 — 2.0 found 
in massive star-forming re gions (Molinari et alj2000l) . Using the 
dust opacities derived by iHildebrandl (1 19831). albeit for grain s 
without ice mantles, as has been done by Beuther et alJ d2002l) . 
leads to masses and column densities that are about a factor of 4 
higher. 

We also list the source multiplicity, in this case all the sources 
above 3 cr found within the 1 1 .4 ' field of view. 
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Table 14. Positions of the LABOCA 870 pm continuum emission peaks. 



Source 




R.A. 


Dec. 


l,b 




Rms 






(J2000) 


J(2000) 






(mJy/beam) 


12326- 


-6245 


12:35:35.0 


-63:02:28.9 


301.14, 


-0.22 


100 


16065- 


-5158 


16:10:19.8 


-52:06:09.7 


330.88, 


-0.37 


50 


16065- 


-5158LABOCA-A 


16:10:23.4 


-52:07:06.1 


330.87, 


-0.38 




16065- 


-5158LABOCA-B 


16:10:07.8 


-52:06:31.4 


330.85, 


-0.35 




16065- 


-5158LABOCA-C 


16:09:55.1 


-52:09:20.8 


330.80, 


-0.36 




16065- 


-5158LABOCA-D 


16:10:44.7 


-52:05:56.2 


330.93, 


-0.41 




16065- 


-5158LABOCA-E 


16:10:47.9 


-52:09:43.0 


330.89, 


-0.46 




16065- 


-5158LABOCA-F 


16:10:21.2 


-52:03:35.9 


330.91, 


-0.34 




16065- 


-5158LABOCA-G 


16:10:15.6 


-52:01:47.4 


330.92, 


-0.31 




16065- 


-5158LABOCA-H 


16:10:13.2 


-52:00:35.5 


330.93, 


-0.29 




16060- 


-5146 


16:09:52.5 


-51:54:55.3 


330.95, 


-0.18 


50 


16060- 


-5146 LABOCA-A 


16:09:49.6 


-51:53:03.7 


330.97, 


-0.15 




16060- 


-5146 LABOCA-B 


16:10:08.7 


-51:58:43.0 


330.94, 


-0.26 




16060- 


-5146 LABOCA-C 


16:10:18.1 


-51:58:43.1 


330.96, 


-0.27 




16060- 


-5146 LABOCA-D 


16:10:17.3 


-51:52:06.2 


331.03, 


-0.19 





Notes. The sub-cores 16065-5158 A-H and 16060-5146 A-D are marked in Fig. HI 



3.2.2. ATCA 

Figure[6]shows the mm continuum emission imaged with ATCA. 
One can see that its distribution agrees well with the peak of 
the LABOCA 870 pm emission. The ATCA continuum emis- 
sion is still unresolved in 16060-5146 and 12326-6245. In 
16060-5146, two cores can be found in the 3 mm emission, 
and one of them, core ATCA-a, is associated with the LABOCA 
peak. The other core, ATCA-b, is too far away from the APEX 
pointing position to be picked up with the APEX beam of about 
19". Therefore, these two cores cannot be responsible for the 
double peaked structure seen in the spectra of 16060-5146. 
Figure [6] also shows the location of the offset positions deter- 
mined for the molecular line observations. While the somewhat 
crude determination of the position via the HCO + peaks meant 
that we may have missed the position of the hot cores by a few 
seconds of arc, we were still close enough to pick up the hot core 
emission in the beams. 

In Fig. [6] overlaid on the data are the interferometric posi- 
tions of associated OH (Caswell et al. 1999, 2001, 2004) and 
CH3OH masers (J. Caswell, private communication), as well 
as the 3 cm and 6 c m peaks of radio continuum emission 
dUrquhart et alj l2007al) . It is remarkable that in 16060-5146, 
only core A (position 16060-5146 ATCA-a) has associated 
masers and radio continuum, suggesting that core 16060-5146 
ATCA-b might be in a very early evolutionary stage. This core 
is neither evident in the LABOCA continuum data nor in the 
Spitzer GLIMPSE MIR data. 

The ATCA data provide a precise position of the core (see 
Table ITSb. The beam de-convolved source sizes, positions, and 
fluxes for the sources were determined from the image plane 
with the imfit task in the MIRIAD software package. 

3.2.3. IR and radio emission 

In addition to the 870 pm continuum observed with LABOCA 
and the 3 mm continuum observed with ATCA (this work), the 
associated fluxes of point sources found in the Spitzer GLIMPSE 
survey in the 3.6, 4.5, 5.6 and 8 pm bands observed with the 
Infrared Array Camera (IR AC) and in the M SX mission at 
8.3, 12.1, 14.7 and 21.3 pm dEgan et al.ll2003h were obtained. 



Photometry on th e Spitzer MIPSGAL 24 and 70 pm images 
dCarev et alj|2006l) could not be performed, because the sources 
are saturated in all cases. All the three regions show extended 
8 pm emission, which is commonly associated with emission of 
photon-dominated regions (PDRs) owing to the IR features of 
polycyclic aromatic hydrocarbons (PAHs). 

In order to find the infrared point sources associated with the 
870 pm dust peaks, we searched the catalogs within a radius of 
9", the LABOCA beam, around the dust peak. When studying 
environments of high-mass star formation, one frequently 
encounters a situation where many infrared sources cluster 
in one single dish beam - in this case the LABOCA beam - 
because massive stars are thought to form exclusively in clusters 
(see Fig . O ■ Because many of the young high-mass proto-stellar 
candidates are very red and bright, they are saturated in the 8 pm 
band and therefore a ppear with null values in the GLIMPSE 
point source catalog dKumar & Gravel [20071) . so when a null 
value at 8 pm was found, we did a manual check on the image 
to distinguish between sources with no detectable emission at 
8 pm and those that are saturated. If the GLIMPSE data was 
saturated at 8 pm, the MSX value was used. 
For 16060-5146 and 12326-6245, both MSX and GLIMPSE 
emission is associated with the hot core traced by the 3 mm 
continuum observed with ATCA. In 16065-5158 however, 
there is neither an MSX nor a GLIMPSE source at the position 
of the hot core, yet both types of sources can be found a few 
seconds of arc offset from the hot core, which might indicate an 
external heating source or a deeply embedded UCH11 region. 

For cm wavelengths, archival 3 cm and 6 cm ATCA data 
published by lUrquhart et al.l d2007al) were obtained for the hot 
cores in 12326-6245 and 16060-5146. In 16065-5158, there is 
an offset of 4.3'.' between the mm position obtained from ATCA 
and the ATCA cm position (see Fig. [6]). The cm emission in 
this source has a negative spectral index of -0.9, hinting at syn- 
chrotron emission. No extragalactic object could be found at this 
position in the NED database. This could be si milar to the situa- 



tion in W3(H20), where a synchrot ron source dReid et alJI 19951) 



can be found in a hot core. Note that lUrquhart et al.T i 2007al) state 
that the spectral indices of individual sources might not be re- 
liable in all cases, because the ATCA is not a scaled array at 
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Table 15. ATCA 3mm continuum parameters 



Source 


Position 


Source Size 


Sy 


S int 




(J2000) 


C) 


(Jy/beam) 


(Jy) 


12326-6245 


12:35:35.06 


-63:02:31.00 


3.3 x 2.0 


1.96 


2.36 


16060-5146 ATCA-a 


16:09:52.57 


-51:54:53.69 


4.6x2.2 


1.24 


1.95 


16060-5146 ATCA-b 


16:09:54.52 


-51:54:53.56 


6.5 x 1.7 


0.78 


1.33 


16065-5158 


16:10:20.00 


-52:06:08.79 


6.4x3.1 


0.15 


0.28 



the observed frequencies (which would provide comparable res- 
olution at the different wavelengths) and the measurements may 
sample different spatial scales. 

Figure [7] shows a zoom-in on the region of the 870 jum peak. 
In 12326-6245, the IR flux is coming from a bright clus- 
ter of objects. One can see emission at 4.5yum (color-coded 
green) surrounding the c luster, which is a sign of shocked gas 
dCyganowski et al.Tl2.Q08b and can point to outflow action. In 
16060-5146, the peak of the 870 fim emission can be associated 
with a single object. There is weak 4.5 ^m emission at the edge 
of the IR source. 16065-5158 has no associated IR source at the 
position of the 870 fim peak. There is however a strong elon- 
gated 4.5 //m emission. This situation is suggestive of a deeply 
embedded young object with a strong outflow. 



4. Derivation of physical parameters 

4.1. LTE modeling 

We observed our sources in lines from CH3CN, CH3OH, and 
H2CO, which are useful tracers of temperature. While classi- 
cally rotational diagrams have been used to determine rotational 
temperature, r rot , and column density, N, of a molecule (see 
iMangum & Woottenldl993h : lQlmi et al.l(ll993h^ we make use of 
an LTE appro ach developed bv | Schilke et al.l ( 1999I) and further 
improved by IComito et al.l (I2005 ). Implemented into the pro- 
gram XCLASS, it allows us to model synthetic spectra that take 
line blends and optical depth effects into account and further- 
more allows us to simulate a double sideband (DSB) spectrum. 
Input parameters for the simulation are source size, temperature, 
column density, line widths, and offset from systemic velocity. 
It is possible to simultaneously model transitions in several fre- 
quency bands and include more than one component, i.e. a core 
and an envelope component and/or several velocity components. 
If there is a hot core and cold envelope component, both com- 
ponents need to be modeled independently though, i.e. the cold 
component is not absorbing the hot component, because we do 
not know the source geometry. Figures. [4] [8] and [9] give an ex- 
ample of a synthetic spectrum for a single molecule overplotted 
on the measured sub-spectra, while Figs.fTOl— IT5lin the appendix 
show the synthetic spectra of all included species. 

The degeneracy between column density and source size 
can be solved for species where both optically thin and thick 
lines are present, as we did for CH3CN and to some de- 
gree CH3OH. For species where only one or two lines were 
available to model, only the beam-averaged column density 
was derived, and for species that we assume to reside in the 
hot c ore, a rotation temperature r, ot =150 K (Schilke et al. 
I2006h was assumed while for species located in the more 
extended envelope 50 K was chosen, unless we had both opti- 
cally thick and thin lines to derive a temperature from the model. 



The lines were identified with the catalogues of the C DMS 
and the JPL dPickett et alj|1998t iMiiller etal]|200U 120051) as a 
reference for the rest frequencies. 

The tables with the detected lines for the three sources can be 
found in the online material (Tables [16] to l28b - here we will 
present the results of the modeling. 

XCLASS produces a synthetic spectrum (model) of all 
species and transitions included in the modeling, which is then 
overlaid on the data. For CH3CN, we performed a y 1 analysis 
on the synthetic model spectrum to determine the best-fit model, 
while for the remaining setups a comparison by eye was done to 
decide on the best fit. This was necessary to account for the high 
complexity when modeling several often blended transitions in 
the double sideband spectra. For sources labeled "b" in Tables 1291 
- |3T] the lines are weak and/or blended with many other lines, 
so that the fit is less accurate than a fit marked "c". 
If both temperature and source size were fixed, the column den- 
sity could be determined to an accuracy of 10%, assuming per- 
fect calibration. The uncertainties introduced when modeling 
both source size and temperature are discussed in the follow- 
ing sections. This modeling is done under the assumption of 
LTE, which may not necessarily be valid for every molecule and 
source, although it is a reasonable choice given the high densi- 
ties. Nonetheless, this method gives us a good and fast overview 
on the chemical composition of the sources, and it can be ap- 
plied to a large sample of sources with considerably less effort 
than individual non-LTE modeling would require. 

The results of the modeling can be found in Tables [29] 
to [3T] Tables [32] to [34] show the model results for those 
species that were modeled in several frequency setups, namely 
CH3OH, CH3CN, H 2 CO and S0 2 . For some species, marked 
as 1 in the Tables, we modeled isotopo logues, using ratios of 
12 C/ I3 C=60, 32 S/ 34 S= 23, 14 N/ 15 N=300 dWilson & Rood! 1 1994 
and 16 O/ 17 O=1500 (IWilson & Rood! 11994 IWouterloot et al I 
120081) . 

4.1.1. CH 3 CN 

The symm etric top molecule C H3CN is an appropriate tempera- 
ture tracer (lArava et al.ll2.005h . The data and the synthetic spec- 
tra are displayed in Fig. [4] Owing to the CH3CN lines up to the 
K=l transition in 16065-5158, we used a compact hot com- 
ponent to model the higher excitation lines and additionally a 
more extended, cooler component to model excess in the K=Q - 
K—2 transitions. Beause this setup also includes some CH' 3 CN 
lines, we had both optically thick and thin transitions to solve 
the N/8 degeneracy and obtain source sizes. While 16065-5158 
and 12326-6245 can both be modeled with a two-component 
fit displaying one hot, compact component and one extended, 
cooler one, it is also possible to obtain reasonable results for 
12326-6245 with only a hot compact component. The two- 
component structure is attributed to a hot dense core inside 
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Table 29. Results of the line modeling for 16065-5158 (-62.2 kms -1 ) 



Species 


Size 




A^b(mol) 




Av 


Offset 


N h /N H2 






(") 


(K) 


(cm" 2 ) 


(cm -2 ) 


(kms" 1 ) 


(kms" 1 ) 






C.H,CN 


a 


150° 


5.0(13) 




3.0 


2.0 


1.7(-10) 


b 


CCH 


a 


50° 


3.2(14) 




6.0 


-1.0 


1.1 (-09) 


c 


CH 3 CCH 


a 


50° 


3.0(15) 




7.0 


0.0 


1.0(-08) 


c 


CH3OCH3 


a 


200 


1.3(15) 




2.5 


3.0 


4.5(-09) 


h 


CH^OCHO-a 


a 


150° 


4.0(14) 




3.0 


0.0 


1.4(-09) 


b 


CO 


a 


50° 


4.0(19) 




5.5 


-0.5 


1.4(-04) 


C J 


CS 


4.0 


90 


3.3(15) 


7.0(16) 


5.0 


0.0 


1.1 (-08) 


C J 


H 2 CS 


a 


50° 


3.0(14) 




6.0 


0.0 


1.0(-09) 


c 


HC 3 N 


0.8 


240 


7.9(13) 


4.0(16) 


7.0 


1.0 


2.7(-10) 


c 




10.0 


40 


4.7(14) 


2.0(15) 


7.0 


1.0 


1.6(-09) 




HCN 


4.0 


150° 


1.4(15) 


3.0(16) 


6.0 


0.0 


4.9(-09) 


C J 


HNCO 


a 


150° 


1.0(14) 




6.0 


1.0 


3.5(-10) 


c 


SO 


8.0 


30 


8.2(15) 


5.0(16) 


7.0 


1.0 


2.9(-08) 


e ,f 



Notes. "Fixed Parameter; *based on weak or partially blended lines only; ''based on one or a few lines only; ^includes vibrationally excited lines; 
•"non-LTE, see text; 'based on isotopologues, see text. The column labeled "offset" shows the offset in velocity relative to the systemic velocity of 
the source. N s are the source- averaged column densities, are the beam-averaged column densities and NJN^ gives beam-averaged abundances. 

Table 30. Results of the line modeling for 16060-5146 (-91 kms" 1 ) 



Species 


Size 


Trot 


iVb(mol) 


jV s (mol) 


Av 


Offset 


N h /N H2 






(") 


(K) 


(cm" 2 ) 


(cm~ 2 ) 


(kms" 1 ) 


(kms" 1 ) 






C 2 H,CN 


10.0 


150° 


2.4(13) 


1.0(14) 


3.0 


0.0 


4.2(-ll) 


b 


CCH 


a 


50° 


5.0(14) 




10.0 


-4.0 


8.8(-10) 


c 


CH 3 CCH 


a 


100 


2.5(15) 




7.5 


-1.50 


4.4(-09) 


c 


CH3OCH3 


a 


150° 


1.0(14) 




2.0 


2.0 


1.8(-10) 


b 


CHiOCHO-a 


a 


150° 


1.0(14) 




2.0 


0.0 


1.8(-10) 


b 


CO 


a 


50° 


6.3(19) 




9.0 


-0.5 


1.1 (-04) 


C J 


CS 


10.0 


55 


1.4(15) 


6.0(15) 


7.0 


-0.5 


2.5(-09) 


C J 


H 2 CS 


a 


50° 


2.6(14) 




10.0 


-0.5 


4.6(-10) 


c 


HC 3 N 


a 


150° 


1.0(13) 




6.0 


0.0 


18(11) 


c 


HNCO 


1.0° 


80 


1.5(14) 


5.0(16) 


6.0 


-2.0 


2.6(-10) 


c 


SO 


4.0 


40 


1.9(16) 


4.0(17) 


6.0 


-0.5 


3.3(-08) 





Notes, indices as in Table 129 1 



a colder, extended envelope. 16060-5146 had to be modeled 
with two hot, compact components with a velocity difference 
of 9 kms" 1 . This effect, which is not obvious in every species 
studied here, suggests a complicated velocity structure for parts 
of the molecular gas in this source. 

We obtained physical sizes assuming near distances for 
16060-5146 and 16065-5158 of 0.03 pc for the hot compo- 
nent in 16065-5158, 12326-6245 and for the two components 
in 16060-5146. The extended components in 16065-5158 and 
12326-6245 are 0.1 pc each. Tables 1321 to [341 show the results 
for the CH3CN modeling. 

4.1.2. CH3OH 

For all three sources, data for the CH30H(7-6) band were ob- 
tained in a setup at 338 GHz, including also the torsionally 
excited v,=l CH 3 OH lines at 337.6 GHz. The CH 3 OH(6-5) 
lines at 290 GHz were observed in all three sources, while for 
12326-6245, the torsionally excited lines at 289 GHz were cov- 
ered in a different frequency setup. Both the 7=6-5 and the 7=7- 
6 series were modeled with XCLASS under the assumption of 
LTE. 



To obtain the temperature for the hot compact component 
in the 7=7-6 band, the optically thin torsionally excited lines 
and the higher excitation lines in the v,-0 bands were used. In 
combination with optically thick lines in the band, the N/9 de- 
generacy could be resolved. The torsionally excited lines at 289 
and 337 GHz also contain a line with higher optical depth, from 
which we could obtain an estimate of the source size for this 
component. 

For details of the model in 12326-6245, see Figs. [8] and [9] 
while plots of the remaining two sources can be found in the ap- 
pendix. The synthetic model spectrum in Fig.[8]overestimates the 
torsionally exited lines, which is because of the added complex- 
ity of modeling the torsionally exited v, = 1 lines simultaneously 
with the hot component of the v, = lines. 

In 16060-5146, the line profiles are much broader owing to 
a blend of two hot compact components and one extended com- 
ponent, and the spectra suffer from a worse baseline than the 
spectra for the other two sources. We could obtain a good fit for 
two hot components, but because of the complexity of the data, 
only a beam-averaged column density could be obtained for the 
cold component. The results of the CH3OH modeling can be 
found in Tables EIHH 
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Table 31. Results of the line modeling for 12326-6245 (-39.3 kms -1 ) 



Species 


Size 


Trot 


A^b(mol) 




Av 


Offset 


N h /N H2 






(") 


(K) 


(cm" 2 ) 


(cm -2 ) 


(kms -1 ) 


(kms -1 ) 






C.H,CN 


a 


150° 


6.0(13) 




3.0 


0.0 


3.1(-10) 


b 


CCH 


a 


50° 


2.8(14) 




4.0 


0.0 


1.4(-09) 


c 


CH 3 CCH 


a 


50° 


1.5(15) 




5.5 


0.0 


7.7(-09) 


c 


CH3OCH3 


a 


70 


1.3(15) 




4.5 


.0 


6.6(-09) 


b 


CH^OCHO-a 


a 


150° 


5.0(14) 




3.0 


0.0 


2.6(-09) 


b 


CO 


a 


50° 


2.8(19) 




4.5 


-0.5 


1.4(-04) 


C J 


CS 


a 


50° 


2.6(15) 




6.0 


0.0 


1.3(-08) 


C J 


DCN 


a 


150° 


3.0(13) 




6.0 


0.0 


1.5(-10) 




H 2 CS 


a 


50° 


2.5(14) 




5.0 


0.0 


1.3(-09) 


c 


HC 3 N 


1.6 


100 


1.6(14) 


2.0(16) 


6.0 


-5.0 


8.2(-10) 


c 


HNCO 


1.6 


75 


2.4(14) 


3.0(16) 


6.0 


0.0 


1.2(-09) 


c 


SO 


10.0 


50° 


7.1(15) 


3.0(16) 


6.0 


0.0 


3.6(-08) 





Notes, indices as in Table 129 1 

Table 32. Properties of 16065-5 158 (-62.2 km s _1 ) as derived from LTE modeling of CH 3 CN, CH 3 OH, H 2 CO and S0 2 . 



Species 


Transition 


Source Size 


Trot 


N s 




Av 


Offset 


N b /N H2 






('■') 


(K) 


(cm -2 ) 


(cm" 2 ) 


(km s" 1 ) 


(kms -1 ) 




CH3CN 


19-18 


0.7 


250 


3.0(15) 


4.5(12) 


4.5 


-1.0 


1.6(-11) 


CH 3 CN 


19-18 


5.0 


30 


1.0(17) 


7.2(15) 


4.5 


-1.0 


2.5(-08) 


CH 3 CN 


16-15 


1.3 


220 


1.4(15) 


7.3(12) 


4.5 


1.0 


2.5(-ll) 


CH 3 CN 


16-15 


5.0 


30 


9.4(15) 


6.7(14) 


4.5 


1.0 


2.3(-09) 


CH3OH 


6-5 


1.8 


140 


2.0(18) 


2.0(16) 


6.0 


0.0 


6.9(-08) 


CH3OH 


6-5 


13.0 


30 


5.0(15) 


1.7(15) 


6.0 


0.0 


6.0(-09) 


CH3OH 


7-6 


1.8 


140 


7.0(17) 


6.9(15) 


6.0 


1.5 


2.4(-08) 


CH3OH 


7-6 


13.0 


30 


5.0(15) 


1.7(15) 


6.0 


1.0 


6.0(-09) 


H 2 CO 


4-3 


2.1 


160 


1.0(16) 


1.3(14) 


6.0 


-0.5 


4.7(-10) 


H 2 CO 


4-3 




50 




4.0(14) 


6.0 


-0.5 


1.4(-09) 


H 2 CO 


6-5 


2.1 


160 


4.8(16) 


6.4(14) 


6.0 


-0.5 


2.2(-09) 


H 2 CO 


6-5 




50 




5.0(13) 


6.0 


-0.5 


1.7(-10) 


so 2 


338 GHz 




30 




3.0(16) 


8.0 


0.0 


1.0(-07) 


so 2 


338 GHz 




250 




1.0(15) 


8.0 


0.0 


3.5(-09) 


so 2 


430 GHz 




30 




3.0(16) 


8.0 


0.0 


1.0(-07) 


so 2 


430 GHz 




250 




1.0(15) 


8.0 


0.0 


3.5(-09) 



Notes. We display in the table the rotational temperature T mt , source-averaged column density N s , beam-averaged column density N b , velocity 
width Av, velocity offset from systemic velocity, and abundances. If two entries per transition are given, they refer to the cold extended envelope 
and the hot compact core. 



4.1.3. H 2 CO 



4.1.4. S-bearing species 



Because it is a slightly as ymmetric rotor, H 2 CO is a g ood tracer 
of kinetic temperatures dMangum & Woottenl Il993[) and it is 
quite ubiquitous in regions of star formation ( Man gum et al] 
ll990HMangum & Woottenlll993l) . 

In this survey, we observed H 2 CO in all three sources in 
the H 2 CO(4-3) transitions at 291 GHz, and in 16065-5158 and 
16060-5146 also in H 2 CO(6-5) lines at 437 GHz. One can ob- 
tain a reasonable fit to the two transitions with a hot compact 
core and a cold extended envelope component. The results of the 
H 2 CO modeling can be found in Tables [3214341 while the spectra 
and synthetic models can be found in Figs. [TO] [12] and[T4]in the 
appendix. For 12326-6245, only the H 2 CO(4-3) transition was 
observed. In 16060-5146, we could not model the hot H 2 CO 
emission with two velocity components, as we did for CH3OH 
and CH3CN. 



We observed the lines from the S-bearing species CS, H 2 CS, 
SO, and S0 2 in our frequency setups. For the H 2 CS(10-9) line, 
only the beam-averaged column density for a fixed temperature 
of 50 K could be determined, while the other species were ob- 
served in several lines. 

We observed the isotopologic C 33 S(9-8), C 33 S(6-5), and 
C 34 S(7-6) transitions in CS, while no C 32 S line was included 
in any of the setups. Assuming isotonic ratios of 32 S/ 3 4 S=23 and 
32 S/ 33 S=127 (IWilson & RooJl994ICh"in et al.ll 19961) . we could 
model CS and SO consistently over the frequency bands. The 
same was not possible for S0 2 . Apart from the S0 2 (5-4) line 
at 351.3 GHz, which has a lower level energy of 19 K above 
ground state, the other S0 2 lines cover higher excitation condi- 
tions, which are between 76 K and 770 K above ground state 
energy. The results of the S0 2 modeling can be found in Tables 
[3214341 while the results for the remaining species are given in 
Tables [29HH] 
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Table 33. Same as in Tabled but for 16060-5146 (-91 km s -1 ). 



Species 


Transition 


Source Size 


Trot 


N s 


N b 


Av 


Offset 


n d /n H2 






('■') 


(K) 


(cm -2 ) 


(cm -2 ) 


(km s _1 ) 


(kms -1 ) 




CH 3 CN 


19-18 


1.1 


170 


1.0(16) 


3.7(13) 


5.5 


2.0 


6.5(-ll) 


CH 3 CN 


19-18 


1.0 


180 


2.0(16) 


6.2(13) 


5.5 


-6.0 


1.1(-10) 


CH 3 CN 


16-15 


1.1 


170 


1.5(16) 


5.6(13) 


5.5 


4.0 


9.8(-ll) 


CH 3 CN 


16-15 


1.0 


180 


3.0(16) 


9.2(13) 


5.5 


-5.0 


1.6(-10) 


CH 3 OH 


7-6 


1.0 


170 


5.0(17) 


1.5(15) 


3.0 


5.0 


2.7(-09) 


CH 3 OH 


7-6 


1.0 


170 


5.0(17) 


1.5(15) 


3.0 


-5.0 


2.7(-09) 


CH 3 OH 


7-6 




30 




4.0(15) 


7.0 


-2.0 


7.0(-09) 


CH 3 OH 


6-5 


1.0 


170 


1.0(17) 


3.1(14) 


3.0 


5.0 


5.4(-10) 


CH 3 OH 


6-5 


1.0 


170 


1.0(17) 


3.1(14) 


3.0 


-5.0 


5.4(-10) 


CH 3 OH 


6-5 




30 




2(15) 


7.0 


-2.0 


3.5(-09) 


H 2 CO 


4-3 


1.3 


170 


5.0(16) 


2.6(14) 


7.0 


0.0 


4.6(-10) 


H 2 CO 


4-3 




70 




3.0(14) 


10.0 


0.0 


5.3(-10) 


H 2 CO 


6-5 


1.3 


170 


5.0(16) 


2.6(14) 


7.0 


0.0 


4.6(-10) 


H 2 CO 


6-5 




70 




6.0(13) 


10.0 


0.0 


l.l(-lO) 


so 2 


338 GHz 




30 




5.0(16) 


10.0 


0.0 


8.8(-08) 


so 2 


430 GHz 




100 




1.0(15) 


8.0 


0.0 


1.8(-09) 



Notes. Because of its line profiles, 16060-5146 was modeled with two hot compact components in some molecules instead of the core/envelope 
model used for the other two sources. 

Table 34. Same as in Tabled but for 12326-6245 (-39.3 km s -1 ). 



Species 


Transition 


Source Size 


Trot 


N, 


N h 


Av 


Offset 


N b /N Hl 






('•') 


(K) 


(cm -2 ) 


(cm -2 ) 


(kms -1 ) 


(kms -1 ) 




CH 3 CN 


19-18/16-15 


1.6 


100 


8(15) 


6.3(13) 


5.4 


0.0 


3.2(-10) 


CH 3 CN 


19-18/16-15 


7.0 


20 


1(15) 


1.3(14) 


5.4 


0.0 


6.7(-10) 


CH 3 OH 


6-5 


1.0 


240 


4(18) 


1.2(16) 


4 


0.5 


6.3(-08) 


CH 3 OH 


6-5 


12.0 


50 


1.5(16) 


4.6(15) 


4 


0.0 


2.4(-08) 


CH 3 OH 


7-6 


1.3 


200 


1(18) 


5.2(15) 


3 


-1.0 


2.7(-08) 


CH 3 OH 


7-6 


12.0 


50 


1(16) 


3.1(15) 


4 


0.0 


1.6(-08) 


H 2 CO 


4-3 




75 




4(14) 


5 


0.0 


2.0(-09) 


so 2 


290 GHz 




35 




2.5(16) 


6 


-0.5 


1.3(-07) 


so 2 


338 GHz 




35 




1.8(16) 


6 


0.5 


9.2(-08) 



4.1.5. Uncertainty estimates 

To estimate the errors of the modeled parameters, we performed 
an x 2 analysis for the CH 3 CN(16-15) band in 12326-6245 
with a fixed source size, varying temperature T, and column 
densities N. This indicated uncertainties of up to 40% in T 
and up to 23% in N within the 3cr confidence limit, which we 
consider as typical for only one frequency setup and no line 
blends from either sideband. It is evident from Tables [3214341 
that in some cases the simultaneous modeling of the hot or cold 
component of a species could not be done consistently in two 
different frequency setups. This discrepancy is smaller than a 
factor of 5 (apart from the cold component of CH 3 CN(19-18) 
in 16065-5158) though, and could have several reasons: 

First, there can be uncertainties in the flux calibration and 
pointing, because the data were taken at different times. In 
addition, in 16060-5158 and 16060-5146 the two setups were 
observed with a positional offset of 8'.'to each other. While 
at these positions the emission is still covered in the beam, 
influences of the source geometry might already play a role. 

Based on the modeling results of the CH 3 CN(16-15) band 
alone, a common line width for the hot and cold component was 
assumed for the CH 3 CN(19-18) transitions. This could be an- 
other potential cause for discrepancies, but the heavy line blend- 



ing of the CH 3 CN(19-18) band does not allow us to derive a 
reliable line width for the hot component. The line widths from 
the XCLASS modeling listed in Tables[32lto[34lare intrinsic line 
widths, and most of the lines are influenced by a high optical 
depth (see also observed average line widths in the results sec- 
tion), which complicates the determination of the line widths. 

4.2. Spectral energy distributions 

To obtain the dust temperatures and luminosities for the three 
sources, the spectral energy distribution (SED; i.e., flux density 
vs. frequency) was analyzed for each source. Estimates of the 
dust temperatures were obtained by fitting graybody plus free- 
free spectra to the SEDs (see Fig. fl6l >. For all three sources 
a power-law dependence of the dust opacity with f3 = 1.8 
was used (see Sec. 13.2. U . based on the dust properties used 
to derive the masses and col umn densities of the three sources 
(lOssenkopf & Henning|1994l) . We also estimated total luminosi- 
ties by integrating the fluxes under the SEDs. Note that the 
GLIMPSE fluxes are not corrected for extinction. 

12326-6245 and 16060-5146 were modeled using three 
different components, while one was used for 16065-5159. The 
data at 3.6 - 12 /^m trace the hot, compact radiation source that 
is already well developed in the infrared, while the data points 
between 24 /mi and 3 mm represent the contribution from the 
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Fig. 2. Maps of the CO(3-2) line for 12326-6245 (top), 
16060-5146 (middle) and 16065-5158 (bottom). The (0,0) po- 
sitions correspond to the 1.2 mm peak positions of Faundez 
(2004) listed in Tab. Q] The velocity and T m b temperature scales 
go from -100 - 20 knX and -5 - 60 K, -150 - knX and 
-5 - 50 K and -120 - knX and -5 - 60 K respectively for 
12326-6245, 16060-5146 and 16065-5158. 



colder, more extended dust envelope. The IRAS points at 60 and 
100 yum are taken from the high-resolution IRAS Galaxy Atlas 
(IGA) data. They represent upper limits to the IRAS flux. The 
cm continuum data, which represent the contribution of free-free 
emission from an associated UCHu region, is modeled as third 
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Fig. 3. C n O(3-2), CO(3-2), CO(4-3) and CO(7-6) transi- 
tions. While C I7 0(3-2) is observed on the peak positions in 
16065-5158 (top), 16060-5146 (middle) and 12326-6245 (bot- 
tom), the spectra of the remaining transitions were produced by 
averaging all the spectra taken at the map positions. The dashed 
line marks the zero level. 



component. In 16065-5158, only the cold, extended envelope 
was modeled, because the radio continuum and infrared data are 
offset from the dust peak. 

For 12326-6245, a dust temperature of 77 K was derived 
with the graybody fit for the extended component. For the hot 
component traced by the Spitzer IRAC points, the blackbody fit 
resulted in a temperature of 390 K. The luminosity of the 77 K 
extended dust component is l.OxlO 5 L Q and the one derived 
from the 390 K compact component 74 L . 
16060-5146 has a dust temperature of 63 K, a temperature 
of the compact component of 500 K and luminosities of 
2.4x10 s L and 490 L Q for the extended and compact compo- 
nents, respectively. 
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12326-6245 



2.941 10 2.942 10 

Rest Frequency (MHz) 



Fig. 4. CH 3 CN(16-15) data at 294 GHz. We show in blue the 
synthetic model spectrum (see Tables |29l- l30l>. Note that the 
spectra for the three sources have been observed with different 
frequency setups, therefore the SO2 line from the image side- 
bands appears at different frequencies. The rest frequencies of 
the CHf CN lines are marked in red. 



For 16065-5158, the dust temperature was derived to be 50 K, 
but in this source, it is hard to constrain the parameters, because 
only the ATCA 3 mm emission and the LABOCA 870 p.m 
emission are associated with it, because IRAS and MSX values 
give but upper limits. The luminosity of the extended component 
is 0.6xl0 5 L o . 

The spectral energy distributions show that the 3 mm sources 
in 12326-6245 and 16060-5146 are dominated by free-free 
emission and by dust radiation for 16065-5158. 



4.3. Spectral types 

To obtain the spectral type of the mm-sources, we esti- 
mated the Lyman con tinuum flux from the cm observations of 
lUrquhart et al.1 (l2007al) . which trace the free-free emission of the 
UCHn region. 

Because that emission is optically thick from both regions, we 
obtained their parameters by fitting the spectral energy distribu- 



301.180° 301.160° 301.140° 301.120° 

Galactic Longitude 
16065-515B 




330.900° 330.650° 

Galactic Longitude 
16060-5146 




331.000° 330.950° 

Galactic Longitude 



Fig. 5. LABOCA 870 fim emission (shown as contours), start- 
ing from 5 <x and continuing in multiples of 1.4 cr (see Table 
[Pil l. The yellow crosses mark the sub-cores in 16060-5146 and 
16065-5158. The LABOCA beam is shown as green circle in 
the lower left of the image. The three-color image shows 3.6 
(blue), 4.5 (green) and 8 //m (red) emission from the Spitzer 
GLIMPSE survey. 
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-52°06'00" I 



-52°06 , 40 l 

16 10 m 22 : 



-5r55'00" I 



-63°02'40' i 




16065-5158 



Fig. 6. ATCA 3 mm continuum data (solid lines) top 
16060-5146 (contour steps -3,3,5,7. ..cr, with cr=0.1 Jy/beam), 
middle 16065-5158 (contour steps -3,3,5,7...<x, with 
cr=8 mJy/beam) and bottom 12326-6245 (contour steps 
-3,3,5,7. ..cr, with cr=0.14 Jy/beam). The synthesized ATCA 
beam is shown in the lower left corner. LABOCA 870 /im 
continuum is shown in gray-scales. Symbols: red: radio contin- 
uum (Urquhart et al. 2007), blue: OH maser positions, green: 
offsets of the APEX line observations, turquoise: IRAS position, 
yellow: CH 3 OH maser (Caswell et al. 1998). 



tion with the th eoretical spectrum of a homogeneous constant 
density plasma dGarav et al.ll2.QQ6l) . For 12326-6245, an emis- 
sion measure of 2.6xl0 9 pc cm" 6 and a Lyman continuum flux 
of 8.2xl0 48 s _1 were obtained, leading to a spectral type between 
Q6.5 and a stel lar luminosity of L — 1.5 x 10 s L , according to 
lPanapaldl973h . 




330.8850° 330.8800° 330.8750° 330.8700° 
Galactic Longitude 

16060-5146 



a -0.1850' 




330.9600° 330.9550° 330.9500° 330.9450° 
Galactic Longitude 

12326-6245 



■H -0.2250' 




301.1450° 301.1400° 301.1350° 301.1300° 301.1250° 
Galactic Longitude 

Fig. 7. Zoom-in of the Spitzer GLIMPSE composite image of 
3.6, 4.5 and 8 fim emission. LABOCA 870 /mi emission is 
shown as contours, starting from 5<x. The small blue circles and 
crosses mark the position of the 6 cm and 3.6 cm continuum 
(Urquhart et al. 2007). 16060-5146: the blue boxes show the lo- 
cation of the ATCA 3 mm continuum emission (marked "a" and 
"b"). 16065-5158: the cross and circle marking the radio con- 
tinuum emission are shown in white. 
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12326-6245 CH30H(7-6), vt=l 
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3.38 1 5 3.38 2 1 5 3.384 1 5 3.38 6 1 5 
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Fig. 8. Spectrum of 12326-6145. The synthetic model spectrum 
of CH30H(7-6) is shown in red. In the upper plot (v t =1), the 
synthetic spectrum of the torsionally excited CH3OH transitions 
(red) is shown. 



12326-6245 CH 3 OH(6-5), vt=l 
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Rest Frequency (MHz) 
12326-6245 CH 3 0H(6-5), vt=0 
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Fig. 9. Spectrum of 12326-6245. The synthetic model spectrum 
of CH30H(6-5) is shown in red. In the upper plot (v t = 1), the 
synthetic spectrum of the torsionally excited CH3OH transitions 
(red) is shown. The two lines marked (mirror) in both plots show 
lines that lie outside the shown frequency range, but were folded 
into the spectrum owing to imperfect filter edges. 



For 16060-5146, the emission measure is 4.0xl0 8 pc cm -6 , the 
Lyman continuum flux 1.5xl0 49 s _1 and the spectral type was 
found to be between 06 and 05.5 with L = 2.5 - 4.0 x 10 5 L Q . 
This kind of analysis works with the assumption, that the Lyman 
continuum luminosity of the source is dominated by the emis- 
sion of one object. 

4.4. Mass estimates 

Table [35] lists the dust masses derived from the 870 //m 
LABOCA observations as well as the virial masses derived 
from the line velocity dispersion of the C 17 Q(3-2) line. M YU 
in M H was derived a ccording to iMacLaren et al.l (1 1988b and 
IMiettinen et alj d2006l) 

Myir = 0.5Av 2 d©, K [M ]. (2) 

The beam-deconvolved source sizes at the 10% flux contour 
were taken as the angular source sizes ©m in k ~ 1.3 is a 
correction factor to correct for the average radial density depen- 
dence. We estimated the dust masses in two ways. M warm gives a 
measure of the masses of the warmer component associated with 
the IRAS source, as traced by the graybody fit of the SED, while 
M ext traces the masses of the extended cooler envelopes. The 
dust m asses were derived using Eq.[3] following Mi ettinen et al.l 
(|2006h. with a dust absorptio n coefficient, of 0.176 m 2 /kg 
( Ossenkopf & Hennindl 19941) (model V) and a dust-to-gas ratio 
^ d = TTio ■ ^v(Td) is the Planck function of a blackbody at dust 
temperature 7d , S v the integrated flux in Jy and d the kinematic 
distance in kpc. 

For the warm component, 7^ as derived from the SED 
and the peak fluxes were used, while the mass of the extended 
dust emission was derived with the integrated intensity within 
the 10% contour and assuming a dust temperature of 20 K 
for 12326-6245 and 30 K for 16065-5158 and 16060-5146, 
respectively. They were assumed to be equal to the gas ki- 
netic temperatures de rived in the CH3CN modeling, following 
IMiettinen et all (120061) . For 16060-5 146, where the CH 3 CN was 
modeled without a cold component, the temperature derived 
from CH3OH was taken instead. 

5. Discussion 

5. 1 . Molecular line data 

The considerable band width of the FFT spectrometers at APEX 
and the huge number of molecular species that we detected in 
each source enabled us to study their chemical composition. To 
facilitate a comparison, the molecular emission was divided into 
two classes, the hot and the cold molecules following the anal- 
ysis described in Sect. 4.1, with rotational temperatures above 
100 K and below 100 K respectively, which corresponds to 
the ice evaporation te mperature for complex organic molecules 
dBisschop et alj|2007l) . 

Several molecules, CH3OH, CH 3 CN, H 2 CO, and S0 2 were 
observed in more than one frequency setup. In these cases it was 
not always possible to model the species with one consistent set 
of temperatures and column densities. The reason for this might 
be either calibration uncertainties between the bands or devia- 
tions from local thermal equilibrium. 
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Table 35. Dust masses derived from LABOCA 870 fim emission. 



Source 




M ext 




**warm 


Rail 


«warm(H 2 ) 


«ext(H 2 ) 




(10 3 M o ) 


(1O 3 M ) 


(1O 3 M ) 


(pc) 


(pc) 


(10 8 cirT 3 ) 


(10 4 crrT 3 ) 


12326-6245 


0.42 


1.6 


3.0 


0.03 


0.6 


6.7 


3.6 


16060-5146 


1.9/5.5 


4.1/12 


13/23 


0.04/0.07 


0.6/1.1 


12/6.9 


6.5/3.7 


16065-5158 


0.51/3.7 


2.0/14 


5.9/17 


0.03/0.08 


0.6/2.1 


8.2/3.1 


1.7/0.63 



Notes. The first entry gives the mass of the hot component traced by the SED, the second of the whole dust envelope traced at 870 fim and the 
third column lists the virial masses derived from C l7 0(3-2) emission. In the last four columns, the radii for the hot compact component (taken 
from the XCLASS line modeling) and the extended dust emission as well as the H2 number densities for both components are listed. If two values 
are listed, they indicate near/far distances used for the calculations. 



In Fig. [T7] the abundances in the hot and cold component are 
shown. In 16060-5146, only CH 3 OH and CH 3 CN were mod- 
eled with two hot components. In the remaining species, it was 
impossible to resolve two hot components, even though large 
line widths often point in this direction. In all three sources, one 
can see a trend of observing higher abundances in the O-bearing 
species than in the N-bearing species in the hot component. 

Chemical differentiation between N- and O-bearing species 
has been observed in several high res olution studies of hot 
molecular cores, among them W3(O H) dWvrowski et alj|1999t) 
and Orion-KL (ISutton et al.1 1 19951) . In the latter, the "hot 
core" is abundant in complex, saturated N-bearing species, 
while the "compact ridge" shows high abundances of com- 
plex Q-bearing molecules. With interferometric observations, 
iLiul (120051) resolve a similar differentiation in the sources E and 
F in G9. 62+0. 19, a high-mass star-forming region with simi- 
lar distance and luminosity to the three sources presented here. 
Because the chemical differentiation is only resolved with high 
spatial resolution observations, one would not expect to find the 
difference in a single dish beam. This leaves the question why 
one sees a higher abundance of O-bearing species in the three 
sources of our sample. One possibility is that the complex O- 
bearing molecules stem from a more extended region and their 
emission is therefore less beam-diluted. 



iRodgers & Charnlevl (1200 ll) argue that a higher abundance 
of O-bearing species might suggest an earlier evolutionary stage 
where the hot core has developed into a rich O-bearing chemistry 
in the first 10 4 yr before evolving into a nitrogen rich state, but 
placing our sources at this young state is inconsistent with the 
high abundance of CH3CN found here. The total fraction of N 
observed is on the order of 3 x 10~ 8 , 3 x 10~ 9 and 5 x 10 -10 
for 16065-5158, 12326-6245 and 16060-5146. This fraction 
was made up of CH 3 CN, HNCO, HC 3 N, and C 2 H 5 CN. Most of 
the nitrogen in hot cores is however locked up in the N-bearing 
species N 2 H + and NH3, which were not covered by any of our 
receiver setups. 

On the other hand, it is conceivable that the complex 
N-bearing species, which according to the above-mentioned 
interferometric studies are expected to be found closest to the 
hot core, have already been destroyed by the radiation of the 
young star, placing the sources at a slightly later stage of the 
hot core phase. For 16060-5146 and 12326-6245, this would 
be consistent with the discovery of cm continuum emission at 
the position of the 870 fim peak. Another origin of the high 
abundance of O-bearing species could be a diff erent composi- 
tion o f the initial grai n mantles as discuss ed by ICharnlev et al.l 
d 19921) and modeled bv lCaselli et~aTI(ll993l) with thermal effects. 



In Fig. [T7] the beam averaged molecular abundances with 
respect to H 2 are compared with those der i ved fo r G327. 3-0.6 
(blue), NGC 6334I (green ISchilke et al.l 120061) and 17233- 
3606 (cyan, iLeurini et al ] 120081) 1 which were also obtained 
with APEX and analyzed in the same way as our data. To 
compare with a well-studied northern hemisphere ho t core, the 
abundances of G34.3+0.15 dMacDonald et al.l 0*9961) were also 
included. Overall, one can see good agreement between our data 
and the other four hot cores for the cold envelope molecules, 
apart from the slight overabundance of SO and S0 2 in our 
sources. 

Comparing the hot species, it is obvious that while they are 
line rich, our sources are weaker in their line emission than 
the other sources. The deficiency of nitrogen-bearing species 
in the hot co mponent is pronounc ed when comparing with the 
other sources. iBeuther et al.l d2009l) do not detect any N-bearing 
species toward two young sources in their high-resolution study 
of hot cores performed with the SMA. They explain the lack of 
N-bearing species by saying that these need more time to evolve 
or warmer gas. Because we are detecting lines from gas with 
T > 200 K, it is not likely that the gas is not warm enough 
for nitrogen chemistry. The association of 16060-5146 and 
12326-6245 with UCHn regions also rules out the hypothesis 
that they are still at an early evolutionary stage. 

A prominent feature of the three sources is the hi gh abun- 
dance of sulfur species. According to I Viti e t al. (2004), the sul- 
fur species are expected to show a time-dependent signature. In 
their models, the authors predict the SO/H 2 CS ratio to be around 
one at earlier times and then lowering, which is not observed in 
any of our sources. 

PVan der Tak et al.l (120031) model the sulfur chemistry in the 
envelopes of young massive stars. The abundances they derive 
for CS and H 2 CS agree with those found in our work. S0 2 and 
SO, however, are observed here at much higher abundances. 
The derived S0 2 abundances compare to t hose of the SQ 2 
"jump " model of the inn er envelopes by Ivan der Tak et al.l 
(l2003l) . lBeuther et al.l (120091) discuss the use of CS, SO and S0 2 
as chemical clocks, with CS as a suitable tracer for the earlier 
stages, while SO and S0 2 better trace the more evolved stages. 
This picture is consistent with our data. 



5.2. Continuum data 

Modeling of the spectral energy distributions of the three 
sources revealed that the 3 mm continuum in 16060-5146 
and 12326-6245 is dominated by free-free emission, pointing 
toward a more evolved stage of the sources. This is consistent 



15 



C. Dedes et al.: three southern HMSF regions 



12326-6245 




1 10 1 10 2 10 3 10 4 10 5 10 6 
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Fig. 16. Spectral energy distributions for the three sources. 
The IRAS fluxes (downward pointing triangles) are taken from 
the high-resolution IRAS Galaxy Atlas (IGA) data. The empty 
squares and triangles show the MSX and Spitzer GLIMPSE 
fluxes respectively, the filled circles the ATCA fluxes and the 
filled triangle and square the bolometer fluxes from SIMBA and 
LABOCA. 
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Fig. 17. Fractional abundances in the three sources for T 
< 100 K (cold component, top plot) and T > 100 K (hot 
component, bottom plot). Red triangles show abundances for 
G34.3+0. 15 (Macdonald et al. 1996), green stars for NGC6334I, 
blue stars for G327.3-0.6 (both Schilke et al. 2006) and turquoise 
stars for 17233-3606 (Leurini et al. 2007). The red lines indicate 
those molecules where the modeling of different transitions re- 
sulted in different column densities (see Tables [34ltol33l. 
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with the infrared data a nd the radio continuum data observed by 
lUrquhart etaP (12007 ah . which also point toward 16065-5158 
as the most embedded and least developed object of the three. 
The luminosities derived from the ionized regions in 
12326-6245 and 16060-5146 differ by about a factor of 
2 f rom the FIR lumi nosities derived from the IRAS fluxes. 
As iGarav etafl (120061) discuss, these differences might result 
because part of the stellar luminosity is absorbed by the dust 
or there might be stars in the region contributing to the FIR 
luminosity, but which are not hot enough to ionize the gas. 
The derived luminosity also indicates that the three sources 
are most likely placed at their near distance, because they 
would otherwise have to harbor an extremely massive star more 
luminous than 04. 

Comparing the dust mass estimates for the three sources, the 
virial and continuum masses agree well for the near distance, 
indicating that the sources are likely bound. To get an estimate 
of the upper limit for the mass that might be concentrated at 
the innermost hot component, we determined the mass using 
the temperature derived from the SED model and the peak 
flux of the 870 jum continuum. These masses are a factor of 
about four lower than the envelope masses for 12326-6245 
and 16065-5158, while they are about a factor of 2 lower for 
16060-5 146, indicating that this source is very centrally peaked. 



Al l three sources have been observed by iFaundez et alj 
(120041) . while 12326-6245 a n d 160 60-51 46 have a ls o bee n 
observed by iMiettinen etall d2006l) and iHill et all d2005h . 
Taking into account the different ways the integrated fluxes 
used, our masses agree well within factors of about 2 with the 
estimates by these authors. 

They are, how ever, more massive th a n the average sourc es in 
the samples o f [ Williams et alj d2004l) . iBeltran et all d2006l) and 
iMueller et al.l ( 2002 ). whose samples contain generally less 
massive and luminous sources. 



6. Conclusions 

We observed the three sources 12326-6245, 16060-5146 and 
16065-5158 in several molecular lines and in the mm and 
sub-mm continuum as a pilot study for a forthcoming line 
and continuum survey of a large sample of 47 high-mass star- 
forming regions. The goal of the pilot study was, first, to define 
and test useful frequency setups to be used with the APEX tele- 
scope for the two surveys and second, to study the three sources 
in more detail, because they all show a hot-core-type molec- 
ular spectrum, which is found in high-mass star-forming regions. 

We chose a setup with HCO + (4-3) and CO(3-2) as our tools 
to pinpoint the location of the dense gas peak as a result of this 
pilot study. These same tools will also provide information about 
the kinematics of the regions from the line shapes of the two 
tracers. The sources were then observed in a setup targeted at 
the CH30H(7-6) series of lines to obtain abundances for the en- 
velope component and a potential hot core component. 
Once LABOCA was operational, the sources were imaged at 
870 //m to study their dust content. 

The sources of the pilot study were observed in five frequency 
setups centered at CH 3 OH(6-5), CH 3 OH(7-6), H 2 CO(4-3), 
H 2 CO(6-5), and CH 3 CN(16-15) bands. The line surveys re- 
sulted in the detection of lin es from 19 different species. 
Following lHatchell et alj d 19981) . the sources are all three clas- 
sified as line-rich sources because of the detection of high- 



excitation CH3OH transitions. Because we also detected com- 
plex organic species, the sources were classified as hot cores. 
We modeled lines from an extended cooler envelope component 
as well as those from a hot compact component. 
All three sources show an overabundance of oxygen-bearing 
species compared to nitrogen-bearing species and a high abun- 
dance of sulfur species. 

16060-5146 and 16065-5158 both show secondary clumps 
in their surroundings. While they are all three located in re- 
gions of active ongoing star formation, the 870 fim peak of 
16065-5158 is not associated with any FIR source. Spectral en- 
ergy distributions including the ATCA 3 mm data were modeled 
for all three sources to derive their dust properties. 
The three sources imaged in 870 fim and 3 mm continuum data 
are massive, luminous hot cores. While 16065-5158 seems to be 
a very young, deeply embedded object at the center of a young 
association, 12326-6245 and 16060-5146 seem more evolved 
and have already developed UCH11 regions. The star-forming ac- 
tivity in the region of 16060-5 146 might have been triggered by 
the large bubble at whose edge it is located. 
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Table 4. Line parameters for CH 3 OH in 16060-5146. 



Frequency 
(MHz) 


Transition 


V 

(km/s) 


5 (v) 
(xlO-'km/s) 


j T mB Av 
(K km/s) 


7 mg 

(K) 


289939.5 
290069 8 


6 0,6^ - 5 QfE 
*J l,6 Aj J 1,5^ 


-90.70 
-90.96 


0.9 
0.4 


17.1 
21.0 


1.5 
2.1 


290110.7 
2901 17 8 


6 ,eA + - 5 ,sA + 

u 5,l Aj ~* 5,0^ 


-91.09 


1.0 


22.1 


2.3 


290183.2 
2901 84 7 


6 4,2-£ ~~ 5 4 ) i£' 
\j 2 5^ — ' 2 4** 


-96.78 


2.5 


8.1 


0.7 


290209.7 
290213.2 


6 3,4^ — 5 3,3-E" 
6 33.E - 5 3,2-E 


-95.49 


10.8 


4.6 


0.5 


290248.8 
290264.2 


6 1 5^ — 5 1 4E 

6 2,4A + - 5 2,iA + 


-91.23 
-89.56 


1.1 
1.7 


16.7 
11.2 


1.4 

0.7 


290307.4 
290307.4 
290307 6 


6 2,5-E - 5 2,4-E 

6 2,5-E - 5 2,4^ 
>J 2,4-*-' 2 3*-^ 


-92.71 
-91.01 


1.6 
0.5 


14.3 
10.0 


1.3 
0.9 


338124.5 

JJOJTt.U 


7 o,7£ - 6 0,6^ 
1 1,7° « 1,6 E 


-91.95 
-92 02 


0.9 
7 


20.4 
29 8 

z.y .0 


1.9 

2 7 


338404.6 
338408 7 


7 6,2^ - 6 g^Zs 

1 0,7^ u 0,6^ 


-92.22 


4.4 


31.3 


2.8 


338456.5 
338475.3 


1 5,2E-6 S .\E 

7 - 6 5,2-E 


-96.36 
-86.17 


2.0 
1.9 


1.9 
1.8 


0.3 
0.3 


338486.3 
338486.3 


7 5 ,2A--6 5 ,iA- 
7 5 , 3 A + -6 5 , 2 A + 


-85.44 


1.4 


3.5 


0.3 


338512.6 
338512.6 
33851 2 9 

_J _JO ~J 1Z..7 


7 4>3 A + - 6 4 ,2A + 
7 4>4 A~ - 6 4, 3 A~ 

' 2 6*^ *^ 2 5*^ 


-87.16 


1.2 


2.1 


0.5 


338540.8 


7 3i5 A + -6 3 , 4 A + 


-80.13 


31.1 


12.4 


0.6 


338559.9 
338583 2 


7 3,5^ - 6 3 , 4 £ 

7 3,4^" — 6 3,3^ 


-91.04 
-92.86 


1.4 
1.2 


18.6 
22.8 


1.2 
1.4 


338615.0 
33861 5 


7 i, 6 £'-6 1,5-E 

' 1,6^ " 1,5 ^ 


-93.69 


2.6 


5.9 


0.6 


338639.9 
338639.9 


7 2 ,5A + -6 2 ,4A + 
7 2 ,5A + -6 2 ,4A + 


-90.94 


0.3 


32.4 


2.8 


338721 6 
338722.9 


7 ^ r F f\ n a1? 

1 2,5- c ' u 2,4- c 

7 2.6-E - 6 2,5-E 


-91 68 

y 1 .vjo 


1.1 


12.6 


1 


431356.4 
CH 3 OH(v,=l) 
337297.4 


9 !, 9 A + - 8 !, 8 A + 
7 lj7 A + -6 lj6 A + 


-91.68 

-95.73 
-86.34 


1.7 

4.7 
3.9 


12.7 

0.6 
1.7 


1.0 

0.2 
0.3 


337605.3 
337610.6 


7 2,5-E - 6 2,4-E 
7 3 , 4 £ — 6 3j3 £ 


-84.84 
-96.86 
-97.77 


3.6 
4.4 
8.4 


0.3 
0.9 
0.5 


0.1 
0.2 
0.2 



Continued on next page, v and 6 (v) are the peak velocity and its error as derived from the Gauss fit. 
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Table 4 - continued from previous page 



Frequency 
(MHz) 


Transition 


V 

(km/s) 


(xlO-'km/s) 


]T mB Av 
(K km/s) 


TmB 

(K) 


337610.6 


7 6,\E - 6 6flE 


Szt Qf\ 




u. / 


\}.L 
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7 2, 5 A + -6 2 , 4 A + 
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-96.48 
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1.8 


0.3 
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7 1,7^-6 i >6 E 
7 ojE - 6 o,6^ 


-98.12 


1.5 


0.6 


0.2 


337655.2 
337655.2 


7 3>4 A- - 6 3l3 A- 
7 3j5 A + -6 3)4 A + 


-84.37 


2.2 


2.0 


0.4 
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Table 5. Line parameters for CH 3 OH in 12326-6245. 



Frequency 
(MHz) 


Transition 


V 

(km/s) 


<5(v) 
(xlO-'km/s) 
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6 4 , 2 A + -5 4 ,iA + 

6 4, 3 A~ - 5 4,2A _ 

6 43E — 5 4 t 2E 


-40.87 


0.7 


23.5 


2.8 


290183.2 

290184 7 


6 4,2-E ~~ 5 

fl n -4" - n A A~ 
y j 2, ~ 2, 4^^ 


-40.77 


5.1 


13.1 


2.4 


290189.5 
290190 5 


6 3 , 4 A + - 5 3 , 3 A + 

u 3,3 rt ~* 3,2^ 


-39.24 


15.3 


12.0 


1.8 


290209.7 
290213 2 


6 3,4-E" — 5 3,3-E 
6 ^ ^3 2-^" 


-39.57 
-39.24 


0.0 
15.3 


7.6 
12.0 


1.2 
1.8 


290248.8 
290264.2 


6 1 — 5 1 4 is 
6 2 ,4A + -5 2 3A + 


-39.08 
-39.07 


1.3 
9.9 


7.4 

13.7 


1.1 
1.7 


290307 6 
290307.6 


A . F S Tip 

6 2 ,4-£' — 5 2j 3-E 


-39 74 


6 5 

VJ. .7 


6 5 


1 


336865 1 

JJUOUJ. J. 

337135.9 
338124.5 

J J O 11 . VJ 


1 2 1 1 1 A~ — 1 2 n n/4~ 

3 ^o-E - 4 2,2 -E 
7 0,7^ - 6 ofiE 

' 1,7^ l,6 i -' 


-39.41 
-39.30 
-39.38 
-39.40 


3.8 
0.3 
0.4 
0.4 


3.5 
16.8 
21.5 
22.1 


0.5 
2.5 
3.2 
3.3 


338404.6 
338408 7 


7 6, 2 £ - 6 

7 mA + -6 n J + 


-39.70 


1.3 


1.5 


0.3 


338442.3 
^nR442 3 


7 6; iA + -6 6;0 A + 

' 6,2^ u 6,1^ 


-39.79 


3.8 


2.1 


0.4 


338456.5 
338475.3 


7 5, 2 £ , -6 5 ,i£' 
7 5>3 £ - 6 5, 2 £ 


-39.61 
-42.32 


1.1 

8.0 


3.5 
15.5 


0.6 
1.8 


338486.3 
338486.3 


7 5;2 A--6 5 ,iA- 
7 5 ;3 A + -6 5 , 2 A + 


-40.86 


8.4 


7.7 


0.9 


338512.6 
338512.6 
338512.9 


7 4 , 3 A + - 6 4 , 2 A + 
7 4j4 A~ - 6 4,3A" 
7 2 , 6 A~ - 6 2,5A" 


-39.78 


1.7 


5.3 


0.8 


338540.8 
338543.2 
338559.9 
338583.2 
338615.0 
338639.9 


7 3 , 5 A + -6 3 , 4 A + 
7 3 , 4 A- - 6 3 , 3 A- 
7 3 ,5^ - 6 3j 4.E 
7 3i 4-E — 6 3>3 -E 
7 i, 6 £-6 i j5 £ 
7 2 ,5A + -6 2 ,4A + 


-40.65 

-39.77 
-39.52 
-38.27 
-39.30 


0.5 

0.8 
0.6 
0.5 
0.7 


24.0 

7.5 
8.7 
27.6 
10.7 


3.4 

1.2 
1.3 
3.4 
1.6 



Continued on next page, v and 6 (v) are the peak velocity and its error as derived from the Gauss fit. 
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Table 5 - continued from previous page 



Frequency 


Transition 


V 


5(v) 


— r 1 

J r, )lB Av 


T m B 


(MHz) 






(km/s) 


(xlO-'km/s) 


(K km/s) 


(K) 


338721.6 


7 2 , 5 £ 


- 6 2,4-E 


-39.91 


0.3 


21.0 


3.0 


338722.9 




- 6 2,5^ 












14 U3 A- 


- 14 0.14A - 




S 7 


1 zt 8 


1 Q 


350103.1 


1 M A + 


-0 ,oA + 


-41.17 


3.4 


1.8 


0.2 


CH 3 OH(v,=l) 














ZOVjjj.U 


6 


- 5 3 3 £ 


/to 1 s 


S 1 

j . 1 


\j.y 


u. 1 


289402.5 


6 2 4^ 


-5 23 £ 


-36.51 


4.8 


2.3 


0.2 


289414.0 


6 24 A + 


-5 23 A + 


-40.00 


19.8 


3.9 


0.4 


289415.0 


6 5 ,2£ 


-5 5 ,iE 










289416.4 


6 4,2-E 


- 5 4J-E 










289420.2 


2,5^ 


- 5 2,4A 


-40.37 


7.6 


1.4 


0.3 


78Qzt7^ S 


6 5 ,i£ 


-5 5 ,o£ 










289427.6 


6 1,6^ 


-5 i, 5 £ 


-38.86 


2.6 


3.6 


0.6 


289428.2 


6 4 , 2 A + 


-5 4 ,iA + 










289428.2 


6 4 , 3 A~ 


-5 41A 










90Q/I9Q 1 


6 3 , 3 A- 


-5 3 ; 2 a- 










289429.1 


6 3 ,4A + 


- 5 3>3 A + 










90Q/19Q o 


6 ofiE 


— 5 0,5^ 










289443.5 


6 2,5-B 


C Z7 

— J 2,4-C 1 


-38.37 


9.7 


2.3 


0.3 


78Q47S fii 


6i, 5 E 


- 5 \fiE 


Q8 
jy.yo 


A 8 


Q 
yj.y 


7 
w.z 


289511.1 


ZT A + 

6 o,eA + 


C A 4- 

- 5 o,5A + 


-40.58 


3.8 


1.4 


0.2 


78Q71 O ^ 

ZO? / 1U.7 


6 t, 5 A- 


-5 i,4A" 






1 Q 
Z. V 


W.J 


337625.7 


7 2 ,sA + 


- 6 2,4A + 


-39.73 


0.9 


1.2 


0.3 


JJ / DJ J.Z 


7 34 A- 


- 6 3 3 A~ 


Q1 


j. 1 


1 n 


O 9 
w.z 


337655.2 


7 3;5 A + 


-6 3;4 A + 










337671.2 


lifiE 




-39.89 


5.2 


1.2 


0.2 


337685.5 


7 4 , 3 A + 


-6 4 , 2 A + 


-39.83 


2.4 


1.3 


0.3 


337685.5 


7 4 , 4 A- 


- 6 4 , 3 A 










337685.5 


7 5,2^ 


-6 5 ,iE 










337707.5 


7 ifiE 


-6 t, 5 £ 


-39.91 


3.4 


1.4 


0.2 


337748.8 


7 o, 7 A + 


-6 , 6 A + 


-39.48 


2.0 


1.3 


0.3 


337969.4 


7 i, 6 A- 


-6 1>5 A- 


-39.63 


1.6 


1.8 


0.4 


350287.7 


15 3,13-E 


- 16 4,i 3 £ 


-38.53 


4.9 


0.8 


0.3 
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Table 6. Line parameters for CH 3 OH in 16065-5158. 



Frequency 
(MHz) 


Transition 


V 

(km/s) 


<5(v) 
(xlO^km/s) 


J T mB Av 
(K km/s) 


TmB 

(K) 


290161 3 


6 4 , 2 A + - 


-5 4 ,iA + 


-61.66 


0.6 


11.4 


1.4 


290183.2 


u 4,2-*^ 


- 5 a i E 


-62.01 


0.4 


14.1 


1.9 


290184 7 


6 o - 

u 2,5^ 


- 5 i,i,4~ 

~* 2,4-** 










290209.7 


6 3,4^ " 


- 5 33-E 


-62.13 


0.3 


15.1 


2.1 


290213.2 


6 3,3ii - 


- 5 3,iE 










290307 4 


6 2 ,s£ - 


-5 2A E 


-62 08 


2 


2 6 


3 


290307.6 


6 2,4^ " 


-5 2 ,iE 










336865 1 


12 U1 A-- 


- 12 o,i2A- 


-61.52 


3.0 


3.2 


0.5 


337135.9 


3 3,o£ - 


-4 2 , 2 £ 


-64.90 


1.6 


15.1 


1.2 


338124.5 


1 nnF - 


u 0,6 ^ 


-64.35 


4.0 


7.7 


0.8 


^^8^44 6 

JJOJTt.U 


7 i, 7 £ - 


-6 l, 6 £ 


-61 89 


1 7 


9 9 


1 .2 


338404.6 


7 tiF - 

' 6,2*-' 


u 6,1^ 


-60.81 


0.7 


5.4 


0.7 


338408 7 


7 o, 7 A + - 


-6 , 6 A + 










338442.3 




- 6 


-61.79 


0.6 


5.8 


0.6 


338442 3 


7 A1y4~ - 


VJ 6,1^ 










338456.5 


7 5 , 2 £ - 


-6 5A E 


-62.08 


0.0 


5.8 


0.9 


338475.3 


' 5,3- c ' 


-6 5 2 E 


-60.00 


10.0 


4.2 


0.7 


338486.3 


7 5;2 A- - 


-6 5 ,iA- 


-61.24 


6.5 


5.1 


0.6 


338486.3 














338512.6 




-6 ^ 


-62.00 


1.0 


12.5 


1.7 


338512.6 


7 4 , 3 A + - 


-6 42 A + 










338512.9 


7 2 , 6 A- - 


- 6 2 , 5 A- 










338540 8 


7 i <4 + - 




-62.20 


0.5 


18.0 


2.4 


338543.2 


7 ^ aA~ - 


-6a 










338559.9 


1 - 

1 3,5^ 


- 6 a /i F 


-61.93 


1.4 


23.1 


2.8 


338583.2 


7 3 , 4 £ - 


-6 3,3£ 


-58.87 


1.5 


1.0 


0.3 


338615.0 


7 i. 6 £ - 


-6 U5 E 


-61.59 


5.1 


2.3 


0.3 


338639.9 




- 6 1 ^/\ + 
u 2,4 /1 


-60.59 


4.4 


3.1 


0.5 


338721.6 




_ 6 -i 


-60.48 


2.6 


2.7 


0.5 


338722.9 


7 2 , 6 £ - 


- 6 










CH 3 OH(v,=l) 














337463.7 


' 6,1** 




-58.92 


2.8 


0.7 


0.2 








-64 48 


3 


5 


2 


337463.7 


7 6, 2 A- - 


-6 6> iA" 










337490.5 


' 6,2*-' 


u b,l L 


-58.05 


2.3 


0.3 


0.1 








-62 63 


5 4 


2 


1 


337519.1 


7 3, 5 £ - 


- 6 3,4E 


-60.14 


3.5 


0.6 


0.1 








-64.34 


1.6 


0.3 


0.2 


337546.0 


7 5,2A" - 


-6 5 ,iA- 


-59.97 


3.0 


1.3 


0.3 








-65.06 


4.5 


0.3 


0.1 


337546.0 


7 5 , 3 A + - 


-6 5;2 A + 










337625.7 


7 2;5 A + - 


-6 2l4 A + 


-58.50 


1.2 


0.7 


0.2 








-61.88 


3.0 


1.0 


0.2 



Continued on next page, v and 6 (v) are the peak velocity and its error as derived from the Gauss fit. 
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Table 6 - continued from previous page 



Frequency 


Transition 


V 


S(y) 


T> 

J T mB Av 


T m B 


(MHz) 




(km/s) 


(xlO-'km/s) 


(K km/s) 


(K) 


J J 1 UJJ. / 


1 2,6^1 O 2,5 A 




Q 6 
y.v 


/t ^ 


yj.j 


337671.2 


7 2fiE - 6 2,5£ 


-61.24 


0.4 


0.6 


0.2 






-65.60 


0.7 


0.4 


0.1 


337685.5 


7 4 , 3 A + - 6 4 ,2A + 


-59.29 


4.0 


0.5 


0.2 




-64.18 


3.7 


0.2 


0.1 


337685.5 


7 4j4 A~ - 6 43A" 










337685.5 


7 5.2S-6 5,1 £ 










337707.5 


7 i i6 £'-6 i, 5 £ 


-60.44 


2.8 


0.8 


0.2 






-63.82 


2.2 


0.2 


0.2 


337748.8 


7„,7A + -6„, 6 A + 


-58.83 


4.1 


0.9 


0.2 






-65.55 


3.4 


0.3 


0.1 


337969.4 


7 i, 6 A--6 i >5 A- 


-60.46 


2.2 


2.5 


0.3 


350287.7 


15 3j i 3 ,E - 16 443.E 


-62.58 


5.2 


0.3 


0.1 
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Table 7. Line parameters for CH 3 CN in 16065-5158. 



Frequency 


Transition 


Velocity v 


S(y) 


J T mB Av 


T,„B 


(MHz) 






(km/s) 


(xlO-'km/s) 


(K km/s) 


(K) 


294025.5 


1 

I67- 


■15 7 


-61.06 


5.7 


3.7 


0.3 


294098.9 


1 a 

166- 


156 


-62.40 


2.9 


7.9 


0.6 


294161.0 


1 A 
ID5- 


1 c 


-59.56 


2.5 


9.0 


0.7 


294251.5 


1 £ 
I63- 


■15 3 


-61.20 


1.0 


11.9 


1.3 


294279.8 


1(32- 


1 c 

"t J2 


-61.27 


3.7 


9.5 


1.1 


294302.4 


1 A 
lOo- 


1 c 

-1 Do 


-58.46 


0.6 


20.1 


1.6 


348911.4 


1 Q 


-109 


-64.58 


1.0 


1.1 


0.2 


349212.3 




1 R, 
-I06 


-62.54 


4.2 


4.2 


0.3 


349426.9 


19 2 - 


■I89 


-61.20 


3.0 


3.5 


0.4 


349453.7 


19 - 


■180 


-61.98 


3.8 


14.8 


0.7 


CH 3 CN(v 8 =l) 














350444.9 


19 - 


-180 


-63.48 


1.5 


1.1 


0.2 


350423.6 


19 2 - 


- 182 


-58.38 


9.6 


0.7 


0.2 


CH' 3 CN 














294151.1 


16,- 


■15, 


-59.85 


1.4 


1.7 


0.3 


294134.2 


16 2 - 


- 15 2 


-61.89 


6.8 


3.4 


0.2 


294106.1 


I63- 


15 3 


-63.40 


1.4 


1.6 


0.4 



Notes, v and 6 (v) are the peak velocity and its error as derived from the Gauss fit. 
Table 8. Line parameters for CH 3 CN in 12326-6245. 



Frequency 


Transition 


Velocity v 


S(y) 


J T mB Av 


T,nB 


(MHz) 






(km/s) 


(xlO-'km/s) 


(K km/s) 


(K) 


294098.9 


16 6 - 


15 6 


-40.11 


4.0 


1.4 


0.2 


294161.0 


I65- 


15 5 


-40.22 


9.4 


1.1 


0.2 


294251.5 


16 3 - 


■15 3 


-39.81 


1.1 


2.3 


0.4 


294279.8 


16 2 - 


■15 2 


-40.04 


0.6 


4.3 


0.7 


294302.4 


I60- 


■15 


-39.88 


3.3 


3.6 


0.6 


349212.3 


19 6 - 


18 6 


-41.06 


0.6 


10.4 


1.4 


349346.3 


19 4 - 


-I84 


-40.43 


3.9 


1.1 


0.2 


349393.3 


19 3 - 


■18 3 


-40.47 


1.7 


1.7 


0.3 


349426.9 


19 2 - 


■18 2 


-40.29 


1.0 


3.8 


0.6 


349453.7 


19 - 


■18 


-40.02 


1.0 


2.9 


0.5 


CH 3 CN(v 8 =l) 














350444.9 


19 - 


18o 


-39.99 


8.9 


0.6 


0.2 



Notes, v and 6 (v) are the peak velocity and its error as derived from the Gauss fit. 
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Table 9. Line parameters for CH 3 CN in 16060-5146. 



Frequency 


Transition 


Velocity v 


S(y) 


J T mB Av 


T m B 


(MHz) 




(km/s) 


(xlO^km/s) 


(K km/s) 


(K) 


294098.9 


16 6 -15 6 


-95.02 


10.1 


2.9 


0.2 






-86.50 


6.4 


0.8 


0.1 


294161.0 


16 5 -15 5 


-96.57 


3.1 


1.3 


0.2 






-86.67 


5.1 


1.8 


0.2 


294211.9 


16 4 -15 4 


-96.28 


6.3 


0.9 


0.2 






-89.23 


13.8 


1.5 


0.2 


294251.5 


16 3 -15 3 


-95.74 


4.6 


3.2 


0.4 






-87.80 


3.8 


2.2 


0.3 


294279.8 


169-15 2 


-92.66 


5.4 


4.7 


0.3 


294302.4 


16o-15 


-88.72 


2.7 


7.3 


0.5 


349212.3 


19 6 -18 6 


-94.13 


5.4 


7.1 


0.3 


349426.9 


W2-I82 


-90.09 


41.5 


1.0 


0.2 


349453.7 


19 ( ,-18o 


-90.03 


3.3 


8.9 


0.6 



Notes, v and 6 (v) are the peak velocity and its error as derived from the Gauss fit. 
Table 10. Line parameters for H 2 CO in 16065-5158. 



Frequency 
(MHz) 


Transition 


Velocity v 
(km/s) 


6(v) 
(xlO-'km/s) 


J T mB Av 
(K km/s) 


TmB 

(K) 


290623.4 
291237.8 


4o,4 - 3 0,3 
42,3 - 32,2 


-62.59 
-61.96 


0.9 
1.3 


27.3 
11.7 


4.3 
1.3 


291380.5 
291384.3 


43,2 - 33,1 
43,1 - 33,0 


-57.40 


5.0 


10.2 


1.3 


436586.5 


62,5 — 52,4 


-61.09 


2.0 


23.3 


2.4 


436751.0 
436751.0 


65,1 - 55,0 

65,2-55,1 


-60.16 


3.2 


13.1 


1.2 


436957.3 
436957.5 


64,3 — 54,2 
64,2 - 54,1 


-61.92 


3.6 


22.1 


1.8 


437199.6 
437236.0 


63,4 - 53,3 
63,3 - 53,2 


-61.92 
-61.51 


3.5 
4.0 


32.6 
36.4 


2.8 
3.1 



Notes, v and 6 (v) are the peak velocity and its error as derived from the Gauss fit. 
Table 11. Line parameters for H 2 CO in 12326-6245. 



Frequency 


Transition 


Velocity v 


6(y) 


J T mB Av 


T m B 


(MHz) 




(km/s) 


(xKT'km/s) 


(K km/s) 


(K) 


290623.4 


4o,4 - 3 0,3 


-39.56 


0.2 


23.0 


3.5 


291237.8 


42,3 - 32,2 


-39.80 


0.3 


12.6 


1.9 


291380.5 


4 3 ,2 - 33,1 


-37.89 


2.3 


14.3 


1.5 


291384.3 


43,1 - 33,0 











Notes, v and 6 (v) are the peak velocity and its error as derived from the Gauss fit. 
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Table 12. Line parameters for H 2 CO in 16060-5146. 



Frequency 
(MHz) 


Transition 


Velocity v 
(km/s) 


( v ) 

(xlO^km/s) 


It a . , 

J ImB&V 

(K km/s) 


T 

I mB 

(K) 


290623.4 
291237.8 


4o,4 - 3 0,3 
42.3 - 32,2 


-92.07 
-92.00 


0.4 
0.7 


26.7 
15.8 


2.7 
1.4 


291380.5 
291384.3 


43,2 - 33,1 
43,1 - 3 3 ,o 


-91.22 


1.2 


14.8 


1.3 


436586.5 


62,5 - 52,4 


-91.65 


2.6 


15.6 


1.2 


436751.0 
436751.0 


65.1 - 55,0 

65.2 -55,1 


-92.69 


10.2 


27.1 


1.6 


436957.3 
436957.5 


64,3 - 54,2 
64,2 - 54,1 


-93.81 


3.0 


12.1 


0.9 


437199.6 
437236.0 


63,4 - 53,3 
63,3 - 53,2 


-93.36 
-93.71 


4.7 
3.0 


15.7 
14.1 


1.3 
1.2 



Notes, v and 6 (v) are the peak velocity and its error as derived from the Gauss fit. 



Table 16. Identified features of CO and HCO + 



v Transition 

(MHz) (J) 

337061.1 3-2 
HCO + 

356734.2 4 , - 3 , 



Table 17. Identified features of CS and HC 3 N 



V 


Transition 


(MHz) 


(J) 


C M S 




337396.5 


7-6 


C 33 S 




291485.9 


6-5 


437165.0 


9-8 


HC 3 N 




291068.4 


32-31 


436468.9 


48-47 


HC 3 N(v7=l) 




291782.3 


32/ -31/ 


337344.7 


37e - 36e 


337825.3 


37/ - 36/ 


437535.9 


48/ -47/ 
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Table 18. Identified features of HCN, DCN and HNCO 



V 


Transition 


(MHz) 




H 13 CN 




431659.8 


5o,o - 4o,o 


HC 15 N 




430235.3 


5o,o - 4o,o 


DCN 




289644.9 


4o,o - 3o,o 


HNCO 




350333.3 


16 1,16 - 15i,is 



Table 19. Identified features of H 2 CO 



v Transition 



(MHz) (] K _, K J 



290623 .4 


4 ,4 


- 3o,3 




291237.8 


4 2 ,3 


- 32,2 




291380.5 


4 3 ,2 


-33,1 


* 


291384.3 


4 3 ,i 


- 33,0 


* 


436586.5 


62,5 


- 52,4 




436751.0 


65,1 


- 55,0 


* 


436751.0 


65,2 


-5 5 ,i 


* 


436957.3 


64,3 


- 54,2 


* 


436957.5 


64,2 


- 54,i 


* 


437199.5 


63,4 


-5 3 ,3 




437236.0 


63,3 


-5 3 ,2 





Notes: (*) Blended with each other 



Table 20. Identified features of CH3OCHO. 



V 


Transition 


Notes 


(MHz) 


(Jjela) 




336889.213 


26o,26£' - 25o,25£ 




336918.095 


266,20^ - 256,19^ 


a) 


338355.771 


27 8 ,i 9 A-26 8 ,i8A 




348909.480 


28 9,20-E - 27g^gE 


b) 


348914.975 


289,2()A - 27g,i9A 




350442.250 


288, 2 q£' - 278, 2 q£' 


c) 



Notes: (a) Blend with CH 3 OH at 350687.8. (b) Blend with CH 3 OH at 337490.5. (c) Blend with CH 3 CN at 350443.3. 
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Table 21. Identified features of CH 3 OH 



V 

(MHz) 


Transition 


Notes 


CH 3 OH(vt=0) 






289939 5 






290069.8 


6 1,6^ - 5 \ t $E 




290110.7 


6 n*A + -5 nsA + 


* 


290117.8 


6 5,1-E - 5 5,0^ 


* 


290138.9 


6 5,2£ - 5 5,i£ 




290145.1 


6 5,1-A" - 5 5i0 A- 


a)* 


290145.1 


6 5 ,2A + -5 5 ,!A + 


* 


290161.3 


6 4 3A" -5 4iA" 


* 


290161.3 


6 4 .2A + -5 4 .iA + 


* 


290162.4 


6 4 3-£" — 5 4 2-E" 


* 


290183.2 


6 4,2 £" - 5 4J-E 


* 


290184.7 


6 o sA~ — 5 o aA~ 


* 


290189.5 


6 3.4A + - 5 3. 3 A + 


* 


290190 5 


a , — 5 -! t4 _ 


* 


290209.7 


6 — 5 33-E 


* 


290213.2 


o 1 1 ft, — i -j or, 


* 


290248.8 


6 1,5^ — 5 1,4-E 




290264.2 


6 2.4A + - 5 2 jA + 




290307 4 


zp 1 ^ 2,4-*-' 


* 


290307.6 


6 2,4-E - 5 2,3-E 


* 


336865.1 


12 1 11 A" - 12 n 19A" 




337135.9 


3 3,0^ - 4 2,2-E 




338124 5 


/ A 1 1*. O A £ 1*. 




338344.6 


7 1,7s -6 i i6 £ 




338404 6 


7 


hi* 


338408.7 


7o.7A + -6o.6A + 


* 


338430.9 


1 \ E — (\ (,<\E 




338442.3 


7 6 ,iA + -6 6 , A + 


* 


338442.3 


7 - 6 * iA~ 


* 


338456.5 


1 * nE — ft s ^ E 




338475.3 


7 5j3 £ - 6 5, 2 £ 




338486.3 


7 5 , 2 A~ -6 5 ,iA" 


* 


338486.3 


7 5,3A + - 6 5 , 2 A + 




338504.1 


7 4 t 4-E — 6 43-E 




338512.6 


7 44 A~ -6 43A" 


* 


338512.6 


7 4 ' 3 A + - 6 4 ^A + 


* 


338512.9 


7 2.6A" - 6 2.5A- 


* 



338530.2 7 4j3 £-6 4 ,2£ 



Continued on next page 
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Table 21. continued. 



Table 21 - continued from previous page 



V 


Transition 

\ } K-,K+> 


Notes 


338540.8 
338543.2 


7 3, 5 A + -6 3A A + 
7 3.4A- - 6 3.3A- 


* 
* 


338559.9 
338583.2 
338615.0 
338639.9 


7 3,5-E — 6 3,4-E 
7 3 t 4-E — 6 3,3-E 

7 1,6-E - 6 h5 E 
7 2 ,sA + - 6 2 .4A + 





338721.6 
338722.9 


1 2,sE -6 2,aE 
1 2fiE - 6 2,5-E 


* 
* 


349107.0 

J JUZo j. j 

350421.6 
350687.7 
430207.1 

^■j 1 jjO.^ 

CH 3 OH(vt=l) 
289355.0 


14 1A3 A~ - 14 

8 „ , 4 8 , „ 4 
2,6^ — i,8 A 

8 i, 7 A-7 2 , 5 A 
4 0,4^ - 3 13.E 
12 2,11-E - 12 1,12^ 

Q „ 4 + 8 „ 4 + 

6 3^-E — 5 3,3-E 

u 2 4-*-' — 2 3 ^ 


d) 


289414.0 
289415.0 
289416.4 


6 2 ,4A + -5 2 ,3A + 

6 5,2^-5 5,l£ 

6 a,iE — 5 4j£ 


* 
* 
* 


289420.2 


6 2,5 A " - 5 2,4A~ 

U 5,l c J 5,0 C 


* 
* 


289427.6 
289428.2 
289428.2 
289429.1 

98QZL9Q 1 
Zo V+ZV. 1 

289429.8 


6 ifiE - 5 1,5-E 
6 4 .2A + -5 4 .iA + 
6 4.3 A - - 5 4 ,iA- 

6 3 ,3A~ - 5 3 ,2A~ 
f\ ^ .4 + S 4 + 

6 - 5 


* 
* 
* 
* 
* 
* 


289475.6 
289511.1 

78Q7 1 ^ 

337297.4 


U 2 5-'-' — 2 4^ 

6 i, 5 £-5 XA E 
6 o.6A + - 5 0.5 A + 

v) ] 5/1 — J 1 4/i 

7 1 ; 7 A + -6 1 > + 




337463.7 


7^,4+ , „ 4 + 

7 6 .2A"-6 6 .iA- 


* 
* 


337490.5 
337519.0 


7 6.2E - 6 e,i£ 

7 3_5-E — 6 3,4-E 




337546.0 


/ 5 9/1 — 05 1/1 

7 5 -?A + - 6 5^A + 


* 
* 


337581.7 


7 4 t 4-E — 6 43-E 


e) 


337605.3 
337610.6 
337610.6 


7 2,5£ - 6 2,4£ 

7 gjis - 6 6,o^ 

7 — 6 3,3-E 


* 
* 
* 


337625.7 
337635.7 


7 2 , 5 A + -6 2 ,4A + 
7 2,eA- - 6 2 , 5 A- 





Continued on next page 
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Table 21. continued. 



Table 21 - continued from previous page 


V 


Transition 


Notes 


(MHz) 


(JiT_,* + ) 




337642.4 


7 ij^-e i i6 £ 


* 




' 0,7 £ - o,6 fc 


* 


337648.2 


7 5j3 £ - 6 5, 2 £ 






/ 3 4A - 33A 


* 


337655.2 


7 3 , 5 A + -6 3 , 4 A + 


* 


337671.2 


7 2,6^ - 6 




337685.5 


7 4. 3 A + - 6 4 . 2 A + 


* 


337685.5 


7 4 , 4 A~ - 6 4,3 A ~ 


* 


337685.5 


7 5,2^ - 6 5,i£ 


* 


337707.5 


7 i >6 £ -6 i <5 E 




337748.8 


7o,7A + -6o,6A + 




337969.4 


7 1.6A--6 1.5A- 




350287.7 


15 3,13^-164,13^ 




13 CH 3 OH 






350103.1 


1 ,,,A + -Oo.oA + 





Notes: (*) Blended with each other, (a) Blend with mirrored line, 
(b) Blend with CH 3 OCHO at 338.396.4. (c) Blend with S0 2 at 338611.8. 
(d) Blend with CH3OH at 338486.3. (e) Blend with 34 SO 337580.1. 
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Table 22. Identified features of S0 2 



V 

(MHz) 


Transition 


Notes 


S0 2 








301896.6 


19? ix - 


- 19, ,q 




338611.8 


20i, 9 - 


-19 28 


a) 


338306.0 


184,4 - 


- 18 3 ,5 


b) 


430193.7 


24 U 3 - 


- 232,22 


c) 


430228.6 


24 ,24 - 


- 23 1,23 


d)* 


430232.3 


212,20 " 


-201,19 


* 


430347.6 


265,21 " 


- 264,22 


e) 


34 so 2 








338320.356 


13 2 ,2- 


-12 u 


b) 



Notes: (a) Blend with CH^OH at 338615. (b) Blend, 
(c) Blend with CH 3 OH at 430207.1 and OCS at 437624.522. 
(d) Blend with HC 15 N at 430235. 
(e) Blend with SO at 430339.5. 



Table 23. Identified features of SO 



V 

(MHz) 


Transition 
(N y ;J,F) 


Notes 


SO 






430339.544 


10io-9 9 


a) 


34 SO 






290562.238 


7 6 -6 5 




337580.147 


8b - 77 


b) 


33 SO 






337198.6 


87 -7 6 





Notes: (a) Blend with 34 S0 2 at 430347.6. (b) Blend with CH 3 OH at 337581.6 



C. Dedes et al.: three southern HMSF regions, Online Material p 15 



Table 24. Identified features of CH3OCH3 and CH 2 CHCN 



V 

(MHz) 


Transition 


Notes 


CH3OCH3 








337420.9 




- 20} ic/t/t 




337421.3 




-20 3 ,i 8 ££ 


* 


337421.8 


2h,i9AE 


-20 3 ,i 8 A£ 


* 


337421.8 


21 2 j 9 £A 


-20 3 ,i 8 £A 


* 


337722.3 


7 44 EE 





* 


337723.0 


7 4A AE 


- 63^^ 




337730.7 


7 4 4 ^4^4 


- 63 3AA 


a) * 


337731.9 


7 4 ' 3 £' J 4 


u3 3^^ 




337732.2 


7 43 EE 


(y^ -^EE 


* 


337770.6 


7 4A EA 


- 6 3 , 4 £A 


b) 


337778.0 


7 \ A EE 


- 63 A EE 


* 


337779.5 


li^AE 


- 6 X4 AE 


* 


337787.2 


li,iAA 


- 63.4AA 


* 


337787.9 


7 43 EE 


- 6?, A EE 


* 


337790.1 


7 4 ' 3 £A 


- &1A EA 


* 


CH 2 CHCN 








337039.730 


362,35 


- 352,34 


c) 



Notes: (*) Blended with each other, (a) Blend with C 2 H 5 OH at 337727.1. (b) bad baseline 
c) Blend with wing of C 17 at 337061.1 
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Table 25. Identified features of CH 3 CN 



V 

(MHz) 
CH 3 CN 
294025.5 
294098.9 
294161.0 
294211.9 
294251.4 
294279.7 


Transition 
I67- 15 7 

166-15 6 

16 5 -15 5 

I64- 15 4 
16 3 -15 3 

162-15 2 


Notes 


294296.7 
294302.4 


I61 - 15i 
I60 - 15 


* 
* 


348911.4 
349125.3 
349212.3 
349286.0 
349346.3 
349393.3 
349426.8 


19 9 -18 9 
19 7 -18 7 
19 6 -18 6 
19 5 -18 5 
19 4 - 18 4 
19 3 -18 3 
19 2 - 18 2 


a) 

b) 
O 


349448.5 
349452.1 


19, - I81 
19o-18 


d) * 

* 


CH 3 CN(v 8 =l) 
350168.1 
350423.6 


19i - 18_i 
19 2 - 18 2 




350444.8 
350449.6 


19 ( ,-18o 
19,-18! 


* 
* 


350465.6 
350507.1 
350552.4 


19 4 - 18 4 
19^ - 18 3 

19 2 -18_2 


e) 


CH* 3 CN 
294066.6 
294106.1 
294134.2 


I64- 15 4 
16 3 -15 3 

162-15 2 




294151.1 
294156.7 


I61 - 15i 
I60 - 15 


* 
f) * 


349220.6 
349254.0 


19 3 -18 3 
19 2 - 18 2 


g) 


349274.1 
349280.7 


19, - I81 
19 -18o 


h) * 

* 



Notes: (*) Blended with each other, (a) Blend with CH 3 OH at 338475.290. (b) Blend with CCH at 349337.7. 
(c) Blend with C 2 H 5 CN at 349392.2 and CCH at 349399.3. 
(d) Blend with C 2 H 5 CN at 349442.9. 
(e) Blend with wing of C 17 at 337061.1. 
(f) Blend with CH 3 CN at 294161.0. (g) Blend with CH 3 OH at 338344.6. 
(h) Blend with CH 3 CN at 349286.0. 
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Table 26. Identified features of CH 3 CCH 



V 


Transition Notes 


(MHz) 


(J) 


290452.2 


17 3 -16 3 


290479.9 


172-16a 


290496.5 


17i-16i 


290502.1 


17o-16 



Table 27. Identified features of CCH 



V 


Transition 


Notes 


(MHz) 


(N,J,F) 





349337.7 4,^ -3,^ a)* 
349339.1 4,|,|-3,|,| * 
349399.3 4,|,|-3,|,| b)* 
349400.7 4,|,|-3,f,f * 

436661.0 5,4,t-4,|,| * 
436661.9 5,^,j-4,J,J * 

436723.1 5,|,|-4,|,| * 
436723.9 5,f,|-4,^| * 

Notes: (*) Blended with each other, (a) Blend with CH 3 CN at 349346.3. 
(b) Blend with CH 3 CN at 349393.3 



Table 28. Identified features of C 2 H 5 CN 



V 


Transition 


Notes 


(MHz) 


Qk.,k + ) 




C 2 H 5 CN(v=0) 






337441.7 


136,17 - 144,10 


* 


337445.9 


374, 33 — 364, 3 2 


* 


349547.0 


399 i3 i - 389, 3 o 


* 


349547.0 


399, 3 o - 38939 


* 


349730.8 


39832 - 38g3i 


* 


349731.3 


398 i3 i - 38830 


* 


349796.0 


39436 - 38435 


a) 


350139.6 


41 1,41 -40i, 40 




350145.1 


41 ,41 -40o,40 




C 2 H 5 CN(v=l) 






349790.2 


30 3 , 8 -29 2 ,7 


a) * 


349791.1 


30 2 ,8 -29 3 ,7 


* 



Notes: (*) Blended with each other, (a) Blend with CH 3 OCH 3 at 349796.6 
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16065-5158 H2C0(4-3) 



16065-5158 H2C0(4-3) 




Rest Frequency (MHz) 

16065-5158 H2C0(4-3) 



Rest Frequency (MHz) 

16065-5158 H2C0(6-5) 




Rest. Frequency (MHz) 

16065-5158 H2C0(6-5) 



Rest Frequency (MHz) 

16065-5158 H2C0(6-5) 




4.372 4.372 10 



Rest Frequency (MHz) 



Rest Frequency (MHz) 



Fig. 10. Molecular emission of 16065-5158 around 290 GHz and 436 GHz. The best-fit synthetic spectrum is overplotted (solid 
line). 
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16065-5158 CH30H(7-6) 



16065-5158 CH30H(7-6) 



■Mi. L , 




3.37 10 3.371 10 3.372 10 3.373 3.373 10 3.374 10 3.375 10 3.376 10 3.377 10 3.378 



Rest Frequency (MHz) 

16065-5158 CH30H(7-6) 



Rest Frequency (MHz) 

16065-5158 CH30H(7-6) 




3.378 10 3.379 10' 



3.38 10 3.381 10 

Rest Frequency (MHz) 



383 3.383 10 



3-384 10 3.385 10 3.386 10 

Rest Frequency (MHz) 



11. Molecular emission of 16065-5158 around 338 GHz. The best-fit synthetic spectrum is overplotted (solid line). 
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16060-5146 H2C0(4-3) 



16060-5146 H2C0(4-3) 
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i 
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l S , 


2.908 10 5 


2.91 10 5 2.912 10 5 




Rest Frequency (MHz) 

16060-5146 H2C0(6-5) 





T„,-M 



2.916 10 
Rest Frequency (MHz) 




4.364 10 4.365 10' 

16060-5146 H2C0(6-5) 



4.366 10 4.367 10 

Rest Frequency (MHz) 



4.368 10 4.369 




Fig. 12. Molecular emission of 16060-5146 around 290 GHz and 436 GHz. The best-fit synthetic spectrum is overplotted (solid 
line). 
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16060-5146 CH30H(7-6) 16060-5146 CH30H(7-6) 




3.378 10 5 3.379 1 5 3.38 1 5 3.381 10 5 3.382 1 5 3.383 3.383 1 5 3.384 1 5 3.385 10 5 3.386 10 5 3.387 1 5 

Rest Frequency (MHz) Rest Frequency (MHz) 



13. Molecular emission of 16060-5146 around 338 GHz. The best-fit synthetic spectrum is overplotted (solid line). 
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12326-6245 H2C0(4-3) 




2.908 10 5 2.91 10 5 2.912 1 5 

Rest Frequency (MHz) 
12326-6245 H2C0(4-3) 



T»bM 





I I 
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- i , S , , i , , , i 



2.914 10 5 2.916 10 5 2.918 10 s 

Rest Frequency (MHz) 



Fig. 14. Molecular emission of 12326-6245 around 290 GHz. The best-fit synthetic spectrum is overplotted (solid line). 
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12326-6245 CH30H(7-6) 12326-6245 CH30H(7-6) 




3.379 10 5 3-38 10 5 3-381 1 5 3-382 1 5 3.383 10 5 3.384 1 5 3.385 10 5 3.386 10 5 3.387 1 5 



Rest Frequency (MHz) Rest Frequency (MHz) 



15. Molecular emission of 12326-6245 around 338 GHz. The best-fit synthetic spectrum is overplotted (solid line). 



